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Preface 


Wafers  of  Ge-Si  (90%-10%)  were  expected  to  have  been  available  at 
RADC/ESO  at  Hanscom  AFB  for  the  fabrication  of  photodetectors  during  the 
summer  of  1979.  Since  this  did  not  come  to  pass,  the  work  of  this  thesis 
focussed  on  the  fabrication  of  Ge  photodiodes  making  use  of  shallow 
junctions  to  increase  the  amount  of  photogenerated  carriers  collected 
by  the  p-n  junction.  This  thesis  describes  the  fabrication  and  testing 
of  ion-implanted  and  diffused  Ge  photodiodes.  It  also  includes  a  quantum 
mechanical  analysis  of  the  energy  band  structure  of  the  proposed  Ge-Si 
alloy  to  predict  the  behavior  of  holes  and  electrons  in  the  crystal  lat¬ 
tice. 

Thanks  are  due  to  Andrew  Yang  for  permission  to  do  the  thesis  work 
at  RADC/ESO.  Special  thanks  go  to  D.  Eirug  Davies,  Joseph  P.  Lorenzo, 
and  Thomas  J.  Ryan  for  their  major  involvement  in  the  project  and  the  ed¬ 
ucation  of  the  author!  Associated  personnel  were  also  very  helpful  during 
the  fabrication,  packaging,  and  testing  procedures  of  each  sample.  I  am 
very  grateful  for  the  assistance  and  direction  provided  by  John  M.  Borky 
in  the  role  of  thesis  advisor.  I  am  forever  indebted  to  my  roommate  and 
friend,  Marla  K.  Doxey,  for  her  typing  skills.  To  my  husband,  Orrin,  I 
dedicate  my  thesis  in  return  for  his  support  and  faith  in  me  during  this 
experience. 


Irene  M.  Nadeau  Mills 
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Abstract 


V>  The  response  of  planar  shallow-junction  Ge  photodiodes  displayed 
peak  performance  at  the  expected  1.46  wavelength.  Responsivity  and 


/  quantum  efficiency  were  measured  at  -vO.9  iikl  fi W  and~80%,  respectively, 
with  leakage  currents  of  ~6  fftC,  typical  for  Ge  photodiodes.  Leakage 
currents  were  notably  decreased  in  diffused  diodes.  The  high  values  for 
y  responsivity  and  quantum  efficiency  could  be  attributed  to  the  shallow 
junction  (y  designed  to  increase  the  amount  of  photogenerated  cur¬ 

rent  whichl could  be  collected  by  the  p-n  junction.  The  reason  for  the 

i 

interest  in  photodetectors  lies  in  fiber-optic  applications.  Minimum 
fiber  attenuation  and  dispersion  occurs  near  1.27  /jm.  Photodetectors  at 
this  wavelength  require  tailoring  of  the  bandgap  energy  of  Ge  by  adding 

i 

about  10%  Si.  A  quantum  mechanical  analysis  of  the  energy  band  structure 
of  the  Ge-Si  alloy  is  included  with  suggestions  to  obtain  energy  gap  and 
effective  mass  values  experimentally.  Ge-Si  photodetectors  fabricated 
as  shallow-junction  photodiodes  should  exhibit  high  photoresponse  and 
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PROPERTIES  OF  ION- IMPLANTED 


AND  DIFFUSED  PHOTODETECTORS  OF 
GERMANIUM  AND  GERMANIUM-SILICON  ALLOYS 
THESIS 


I  Introduction 


Background 

Within  the  last  decade,  communication  using  light  signals  transmitted 
over  optical  fibers  has  become  a  feasible  alternative  to  the  transmission 
of  electrical  signals  over  wire.  Optical  fiber  systems  are  attractive 
for  Air  Force  applications  because  of  their  advantages  over  present  sys¬ 
tems:  elimination  of  crosstalk,  reduced  size  and  weight,  lower  cost, 
low  losses  and  high  bandwidths,  plus  freedom  from  ground  loop  pickups 
and  electromagnetic  interference  (EMI)  (Ref  38:5-6,10).  The  development 
of  light  sources,  particularly  for  use  in  fiber  communication  systems, 
has  spurred  the  development  of  photodetectors.  The  light  sources  employ¬ 
ing  III-V  semiconductor  materials  include  ternary  and  quaternary  alloys 
whose  variable  composition  makes  it  possible  to  fabricate  sources  for  a 
wide  span  of  wavelengths  (Ref  12:110).  Light  sources,  such  as  lasers  and 
light-emitting  diodes,  convert  electrical  signals  into  optical  radiation 
In  wavelengths  that  range  from  the  far  infrared  (1-100  p m)  region  through 
the  visible  (0.45-0.70  pm)  and  the  ultraviolet  (0.01-0.1  pm)  regions. 

The  optical  radiation  transmitted  through  an  optical  fiber  is  converted 
back  to  electrical  signals  for  signal  processing  in  the  receiver  by  a 
photodetector. 
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Fig  1.  Spectral  attenuation  curve  of  optical  fiber  with  8%  GeO„  -  0.7% 
doped  Si02  core(Ref  58) . 


In  looking  at  desirable  optimization  of  the  source-channel-detector 
system  (Ref  3;  28;  56),  certain  fiber  characteristics  dictate  preferred 
wavelengths  of  operation,  that  is,  the  wavelengths  at  which  transmission 
through  the  fiber  channel  will  result  in  minimum  losses.  Refer  to  the 
spectral  attenuation  curve  of  a  typical  fiber  in  Figure  1.  In  consider¬ 
ing  three  wavelengths  of  operation:  ~0.85  am,  ~1.0  am,  and  —  1 . 27  am, 
some  interesting  points  have  been  observed  (Ref  12).  At  approximately 
0.85  am,  the  popular  wavelength  at  present,  GaAlAs  sources  are  well  match¬ 
ed  to  Si  detectors.  At  ~1.0  am,  silicon  detectors  experience  a  drop  in 
quantum  efficiency  when  matched  to  GalnAsP  LED's,  injection  lasers,  and 
Nd:YAG  crystals;  they  also  require  a  large  applied  bias  sufficient  to 
obtain  the  needed  100-200  am  depletion  layer  width,  increasing  the  leakage 


current,  and  hence  noise.  The  interest  in  the  1.0  n m  wavelength  is  due 
to  reduced  fiber  attenuation  and  dispersion,  that  is,  improved  fiber 
characteristics  that  help  compensate  for  the  shortcomings  of  the  Si  de¬ 
tector  at  1.0  (!m.  Absorptionless  transmission  is  possible  and  material 
dispersion  in  the  fiber  is  minimized  for  wavelengths  near  1.27  //m 
(Ref  12;  58).  However,  because  silicon  detectors  are  not  sufficiently 
sensitive  at  1.2-1. 3  /m,  interest  has  shifted  to  other  III-V  compounds 
to  fabricate  detectors  which  have  greater  quantum  efficiency.  In  the 
near  1.3  fim  range,  III-V  compound  sources  are  presently  available.  Good 
detectors  for  operation  at  1.27  are  the  next  hurdle.  It  is  well  known 
that  the  wavelength  at  which  a  semiconductor  material  will  detect  is  a 
function  of  its  energy  gap.  This  implies  that  a  semiconductor  crystal 
must  be  grown  whose  energy  gap  brings  about  photosensitivity  at  a  wave¬ 
length  of  ~  1 . 27  #i m. 

Status  of  Germanium  Diodes 

Germanium  has  the  ideal  bandgap  (0.67  eV)  for  detecting  wavelengths 
in  the  range  of  1.4-1. 6  */m  as  shown  in  the  spectral  response  curves  of 
Si  and  Ge  in  Figure  2  (Ref  34).  It  does  detect  wavelengths  between  1.2- 
1.3  /jm  but  an  energy  gap  of  about  0.95  eV  would  be  ideal  here.  One  pos¬ 
sible  approach  is  the  use  of  III-V  compound  alloys  to  "tailor"  the  energy 
gap.  The  Ge-Si  system  shown  in  Figure  3  suggests  that  the  energy  gap  of 
Ge  can  be  tailored  by  combination  with  Si  to  meet  the  requirements  of 
detection  at  ~1.27  fim.  It  is  proposed  that  a  Ge-Si  alloy  would  produce 
the  correct  bandgap  for  detecting  wavelengths  which  the  fiber  transmits 
best.  Energy  bandgap  is  drawn  as  a  function  of  the  percentage  of  Si  in 
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Fig  2.  The  spectral  response  Curves  of  S.i  and  Ge  photodiodes  (Ref  34) . 

Ge  and  is  scaled  on  the  left;  the  lattice  constant  scale  is  on  the  right. 
The  graph  shows  that  the  bandgap  does  not  increase  linearly  as  the  per¬ 
centage  of  Si  in  Ge  increases.  In  fact,  the  function  has  a  break  in  it 
showing  a  steep  slope  when  0-15%  Si  is  added  to  Ge  and  a  less  steep  slope 
elsewhere.  The  steeper  slope  implies  that  a  larger  variation  in  bandgap 
is  possible  with  a  small  amount  of  Si  in  Ge.  However,  it  should  be  sta¬ 
ted  that  the  mismatch  in  the  lattice  constants  of  Ge  and  ~10%  Si  begins 
to  cause  a  breakdown  in  the  crystal  lattice.  But  according  to  the  curve, 
a  10%  Si  -  90%  Ge  alloy  suggests  a  bandgap  of  ~0.85  eV,  based  on  the  re¬ 
lationship  between  energy  gap  and  wavelength  of  detection,  which  approach¬ 
es  the  desired  bandgap  for  detection  at  a  wavelength  near  1.27  flm. 
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Fig  3.  The  Ge-Si  system  (Ref  5;  27). 


Problem  and  Scope 

One  of  the  major  problems  in  the  proposal  to  fabricate  semiconductor 
photodetectors  from  a  Ge-Si  crystal  alloy  is  the  growth  of  the  crystal 
alloy  itself.  The  change  in  interest  from  Ge  to  Si  in  the  early  days  of 
semiconductor  technology  has  almost  made  Ge  a  forgotten  material.  Much 
of  the  research  and  study  of  Ge  and  Ge-Si  alloys  was  dropped  when  the 
convenience  and  versatility  of  the  Si-SiC^  system  was  realized.  It  is 
necessary  to  study  the  early  literature  on  Ge,  especially  Ge-Si  alloys, 
and  attempt  to  grow  the  alloy  using  the  modern  and  well  defined  techniques 
for  growing  crystals.  In  the  meantime,  while  progress  is  being  made  by 
crystallographers,  the  problem  of  fabricating  devices  from  Ge-Si  can  be 
addressed . 

The  first  step  toward  fabricating  devices  from  the  alloy  is  to  work 
with  the  basic  substrate  material,  ie.  Ge.  That  is  the  main  subject  of 
this  thesis.  In  other  words,  since  the  crystal  alloy  of  interest  is 
primarily  Ge,  an  indication  of  the  behavior  of  devices  made  in  Ge-Si  can 
be  obtained  from  those  made  in  pure  Ge.  The  work  involved  in  this  thesis 
is  limited  to  Ge  photodiodes  with  very  shallow  junctions.  The  goal  is 
to  make  very  shallow  junction  Ge  photodiodes  by  well  known  fabrication 
techniques  such  as  diffusion  and  ion  implantation  which  are  comparable 
in  performance  to  Si  photodiodes,  each  in  its  own  preferred  wavelength 
range  of  operation.  A  more  immediate  goal  is  to  verify  the  performance 
of  the  Ge  photodiodes  with  known  results  such  as  peak  photoresponse  at 
~1.46//m  expected  of  Ge  and  responsivities  comparable  to  other  Ge  photo¬ 
diode  structures.  Once  techniques  are  developed  for  making  Ge  photodiodes 
which  work  reasonably  well,  the  next  step  will  be  to  tailor  the  bandgap 
of  the  Ge  substrate  material  with  Si,  thereby  changing  the  wavelength  of 
operation. 
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General  Approach 


In  the  interest  of  obtaining  photodetector  performance  at  1.27 //m, 
the  approach  of  concern  in  this  thesis  begins  with  Ge  photodiodes  whose 
peak  response  occurs  at  ~1.46/im.  The  Ge  diodes  which  are  fabricated 
must  have  increased  quantum  efficiency  and  decreased  leakage  currents 
from  the  large  values  which  have  been  characteristic  of  Ge  devices.  The 
distinguishing  feature  between  other  Ge  photodiodes  and  those  pursued  in 
this  work  is  the  shallow  junction,  —  0.25-0.75  /Zm  deep,  discussed  further 
in  Section  II.  The  goal  of  diffused  diodes  is  uniform  doping  across  the 
active  area  of  the  diode.  The  diffused  diodes  are  mesa  etched  to  reduce 
the  surface  leakage  currents  and  define  the  active  area  of  the  diode. 

The  ion  implanted  diodes  promise  better  control  for  obtaining  shallow 
junction  depths  and  doping  concentrations.  Implanted  diodes  discussed 
in  this  thesis  have  been  thermally  annealed,  although  there  is  interest 
in  laser  and  electron  beam  annealing. 

Each  completed  sample  photodiode  was  tested  in  two  areas  of  perform¬ 
ance.  First,  each  sample  was  tested  as  a  diode  to  verify  expected  diode 
action  and  for  low  leakage  currents.  Second,  each  diode  is  tested  for 
its  optical  properties.  The  photoresponse  characteristics  determine  the 
wavelength  of  peak  response  and  the  quantum  efficiency  of  the  photodiodes 
which  can  be  compared  to  previously  reported  Ge  photodiode  photoresponse 
characteristics  as  listed  in  Table  I. 

Work  done  in  parallel  with  the  fabrication  of  Ge  photodiodes  includes 
the  crystal  growth  of  the  Ge-Si  alloy  that  will  supply  the  material  need¬ 
ed  for  devices  optimized  for  a  wavelength  of  operation  of  1.27  //m.  The 
first  step  is  to  be  able  to  grow  an  epitaxial  layer  of  Ge  on  top  of  a 
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known  Ge  sample.  The  ratio  of  Si  to  Ge  is  then  increased  slowly  to  en¬ 
courage  formation  of  single  crystal  Ge-Si  alloy.  Once  such  an  alloy 
exists,  fabrication  of  the  alloy  photodiodes  using  the  techniques  devel¬ 
oped  for  pure  Ge  substrate  photodiodes  can  be  attempted. 

Sequence  of  Presentation 

With  the  forgoing  discussion  as  an  explanation  of  the  interest  in 
Ge  and  Ge-Si  photodiodes,  Chapter  II  provides  a  background  to  the  general 
area  of  photodiodes  and  how  they  operate.  The  discussion  includes  simple 
photodetector  theory  to  provide  motivation  for  shallow  junctions.  This 
is  followed  by  a  discussion  of  the  diode  characteristics  that  result  from 
the  fabrication  techniques  used  to  form  these  shallow  junctions.  Chapter 
III  gives  the  theoretical  background  for  tailoring  the  bandgap  of  semi¬ 
conductor  materials  and  the  effects  of  this  tailoring  on  wavelength  of 
response,  specifically  of  the  Ge-Si  system  to  obtain  peak  response  at 
1.27  //m.  Chapter  IV  provides  a  history  of  the  fabrication  of  each  sam¬ 
ple  tested  and  Chapter  V  outlines  the  experimental  procedure  and  equip¬ 
ment  used. 

The  data  reduction  procedure  and  results  are  given  in  Chapter  VI. 
Limits  of  the  validity  of  the  data  are  also  discussed  here.  An  analysis 
of  the  results  is  given  in  Chapter  VII.  It  covers  photoresponse  and 
diode  characteristics  of  the  implanted  and  diffused  diodes.  Conclusions 
and  recommendations  are  presented  in  Chapter  VIII. 
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II  Background 


Introduction 

Coherent  and  incoherent  light  sources  such  as  lasers  and  LED's  pro¬ 
duce  wavelengths  of  optical  radiation  that  extend  from  the  far  infrared 
region  through  the  visible  to  the  ultraviolet  region.  For  this  reason, 
photodetectors  are  being  developed  with  high  sensitivity  to  weak  light 
signals  and  response  to  light  intensity  modulation  at  frequencies  extend¬ 
ing  into  the  microwave  region  (Ref  36;  37) .  The  development  of  various 
light  sources  has  subsequently  spurred  the  development  of  photodetectors. 
It  is  possible  to  design  optical  detectors  and  optical  receivers  for 
performance  at  various  wavelengths  of  light  and  combinations  of  speed  of 
response. 

An  understanding  of  the  photoelectric  effect  provides  insight  into 
how  photodetectors  operate  at  certain  wavelengths.  The  speed  of  response 
of  photodetectors,  i.e.  the  bandwidth  of  modulating  signals  which  can  be 
processed,  depends  on  other  processes  associated  with  semiconductor  de¬ 
vices. 

The  Photoelectric  Effect 

The  generation  of  carriers  in  a  semiconductor  by  incident  light  is 
based  on  the  photoelectric  effect  (Ref  15:31-38).  According  to  this 
phenomenon,  a  photon  with  energy 

E  =  YiV  (1) 
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where 

E  =  energy,  eV 

3  A 

h  =  Planck's  constant  6.625  x  10  joule-sec 
V  =  frequency,  Hz 

which  is  inherent  in  the  light  incident  on  the  photodetector  is  absorbed 
by  the  semiconductor  material  and  forces  the  transition  of  a  carrier  from 
a  quantum  state  of  low  mobility  to  a  quantum  state  of  high  mobility.  The 
conductivity  of  a  material  is  a  function  of  carrier  mobility  and  the  num¬ 
ber  of  carriers  as  expressed  by 

a  =  q  (  p)  (2) 

where 

O  =  conductivity,  1/42  . 

-19 

q  =  charge  of  the  carrier,  1.602  x  10  coul 

2 

=  mobility  of  electrons,  cm  /V-sec 

_3 

n  =  density  of  electrons,  cm 

2 

/ip  =  mobility  of  holes,  cm  /V-sec 

_3 

p  =  density  of  holes,  cm 

When  incident  light  falls  on  the  surface  of  the  photodetector,  carriers 
are  generated  either  by  intrinsic  or  extrinsic  transitions  (Ref  53:628- 
629).  In  the  case  of  intrinsic  transitions,  the  increase  in  conductivity 
under  illumination  is  due  to  the  increase  in  the  number  of  carriers. 

The  wavelength  at  which  detection  occurs  is  derived  from  the  photon  energy 
and  the  energy  bandgap  of  the  semiconductor  material.  When  the  amount 
of  photon  energy  is  just  sufficient  to  cause  the  transition  of  a  carrier 
from  its  bound  state  in  the  lattice  to  a  free  state,  that  photon  energy 
is  equal  to  the  energy  bandgap  of  the  material. 
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Wavelength  is  related  to  frequency  through  the  velocity  of  light, 
c  =  \v  .  Substituting  for  v  into  Equation  (1): 

Eg  =  hv  =  ~  (3) 

or 


where 


X 


he 

E 

g 


1.24 

E  (eV) 
6 


jiin 


(4) 


X  -  wavelength,  //m 
h  =  Planck's  constant 

g 

c  =  velocity  of  light,  3  x  10  m/s 

E  =  energy  bandgap,  eV. 

g 

The  important  result  of  this  simple  substitution  is  that  the  wavelength 
which  a  semiconductor  material  detects  is  related  to  the  energy  gap,  E^, 
of  that  semiconductor  material.  Therefore,  to  detect  optical  radiation 
at  a  particular  wavelength  or  range  of  wavelengths,  one  must  consider 
the  material  the  photodetector  is  made  of. 

It  is  also  possible  to  have  energy  states  within  the  forbidden  energy 
gap  and  for  transitions  to  occur  between  these  states  and  the  higher 
mobility  conduction  states.  This  process  is  known  as  extrinsic  photo¬ 
excitation.  By  means  of  doping  with  impurities  whose  allowed  energy  states 
fall  in  the  forbidden  energy  gap  of  the  parent  material,  the  device  can 
be  made  to  detect  optical  wavelengths  longer  than  those  predicted  by 
Equation  (4).  Conductivity,  as  a  result  of  photoexcitation,  takes  place 
through  the  absorption  of  photons  of  energy  greater  than  or  equal  to  the 
energy  separation  of  the  conduction  and  valence  band  levels. 
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Shallow  Junctions 


The  current  that  can  be  generated  in  a  photodiode  is  dependent  on 
the  amount  of  light  absorbed  by  the  material.  The  absorption  coefficient, 
a  (cm  ^) ,  is  also  a  strong  function  of  wavelength  and  semiconductor 
material  as  shown  in  Figure  4  for  Si  and  Ge  at  300°K  (Ref  14:1151;  53:661). 
Germanium  can  absorb  light  in  a  range  of  wavelengths  between  0.65  and 
1.62  fim.  The  absorption  coefficient  and  wavelength  are  factors  that 
determine  the  intensity  of  light  as  it  penetrates  the  semiconductor, 
according  to  the  relationship 

I  =  IQ  exp(-ax)  (5) 

where 

I  =  light  intensity  at  a  depth,  x,  into  the  semiconductor  from 
the  surface 

Ip  =  total  incident  intensity 
Cl  =  absorption  coefficient,  Ot  (  X)  cm  * 
x  =  distance,  cm. 

The  intensity  falls  off  exponentially  with  distance  into  the  material. 

A  plot  of  intensity  versus  distance  is  shown  in  Figure  5  for  wavelengths 
0.65,  1.0,  1.27,  and  1.5 flm,  at  each  of  which  there  exists  some  interest 
for  optical  fiber  systems.  The  absorption  coefficients  for  these  wave¬ 
lengths  are  10^,  6.6  x  lo\  8.0  x  10^,  and  5.0  x  10^  cm  *,  respectively, 
as  obtained  from  Figure  4.  Certain  deductions  can  be  made  by  observing 
the  curves  in  Figure  5.  A  discussion  of  this  follows. 

It  is  possible  to  determine  the  depth  and  width  of  the  depletion 
region  needed  for  absorption  of  specific  wavelengths.  Only  carriers 
generated  within  the  depletion  region  and  within  one  diffusion  length 
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Fig  5.  Exponential  decay  of  light  as  it  is  absorbed  into  a  material. 

of  the  depletion  region  can  be  swept  out  by  the  field  in  a  p-n  junction 
and  be  collected,  thus  contributing  to  the  detected  signal  current.  For 
example,  very  shallow  junctions  are  essential  for  collecting  light  at 
short  wavelengths  and  little  applied  bias  is  needed  to  sufficiently  widen 
the  depletion  region. 

It  Is  obvious  why  Si  is  chosen  over  Ge  as  the  photodetector  sub¬ 
strate  for  devices  operating  at  0.65  //m.  Si  has  an  absorption  coefficient 

3 

about  two  orders  of  magnitude  lower  than  Ge  at  this  wavelength  (3  x  10 
cm  *  for  Si  versus  10"*  cm  *  for  Ge).  The  majority  of  the  light  which 
managed  to  penetrate  the  Ge  substrate  would  be  consumed  by  window  losses 


* 
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which  will  be  explained  later.  The  absorption  coefficient  at  1.0  tt m 
was  chosen  for  comparison  because  silicon  "cuts  off"  at  that  point.  The 
1.27  nm  value  of  absorption  coefficient  was  chosen  for  consideration 
because  of  the  significance  it  has  in  fiber  characteristics  as  a  region 
of  minimum  dispersion  and  a  fiber  attenuation. 

Figure  6  shows  a  plot  of  the  depletion-layer  width,  V(flm)  as  a 

_3 

fjnction  of  substrate  doping,  N  (cm  )  and  applied  reverse  bias.  It  is 

D 

obtained  by  assuming  the  profile  of  a  shallow  junction  approximates  a 
one-sided  abrupt  junction  (Ref  53:89)  and  using  the  equation 


where 

W  ■  total  depletion  width,  fi m 

«  ■  semiconductor  permittivity,  farad/cm  (16  for  Ge) 

Vbi  =  built-in  potential,  volts 

Vbias  =  aPPH-e<*  bias  voltage,  volts. 

The  wafers  used  to  make  the  devices  studied  in  this  thesis  were  3  Q-cm 

15  -3 

p-type  Ge  wafers,  for  which  Nc  ~10  cm  .  Using  Figure  6  in  conjunction 

D 

with  Figure  5,  with  the  assumptions  of  a  ~0.75  #/m  junction  depth  and  at 
least  a  3  one-sided  depletion-layer  width,  shows  that  the  device 
could  capture  about  50%  of  the  incident  light.  At  1.5  n m,  the  absorption 
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BACKGROUND  CONCENTRATION,  NB  (cm“^) 

a  function  of  substrate  doping  and  applied  reverse  bias. 


coefficient  for  Ge  drops  off  rapidly.  Peak  response  has  been  observed 
to  occur  near  1.5  /im.  The  corresponding  line  of  Figure  5  shows  that 
1.5  /im  radiation  will  penetrate  deep  into  the  device,  far  beyond  a  single 
diffusion  length  from  the  depletion-layer  width.  The  optically-generated 
carriers  would  be  recombined  in  the  bulk  and  never  have  a  chance  to  be 
collected  at  the  p-n  junction.  To  reemphasize  the  focus  of  this  thesis, 
at  1.27  pm,  the  efficient  absorption  of  radiation  would  require  a  shallow 
junction  depth  in  the  range  0.25-0.75  pm  and  a  depletion-layer  width  of 
~  5.5  /jm. 

Photodetector  Losses 

Under  reverse  bias,  the  depletion  region  of  a  p-n  photodiode  widens 
into  the  lightly  doped  side  because  of  the  one-sided  abrupt  junction 
assumption  as  shown  in  Figure  7.  Also  shown,  are  the  various  losses  that 
are  incurred  by  the  viewpoint  of  a  photodetector.  The  window  loss  is 
denoted  by  light  wave  (1).  The  n-type  As  doping  is  the  window  through 
which  the  radiation  must  penetrate.  Absorption  in  the  oxide  or  oxide 
equivalent  is  denoted  by  (2) .  To  minimize  absorption  and  reflection 
losses  an  antireflection  (AR)  coating  which  exhibits  the  quarter-wave 
plate  effect  may  be  needed.  The  third  type  of  loss,  (3),  shows  the  path 
of  radiation  as  it  penetrates  deep  into  the  p-type  Ge.  The  e-h  pairs 
that  result  from  carrier  generation  have  a  greater  probability  of  being 
recombined  before  drifting  back  to  the  junction  where  they  can  be  col¬ 
lected.  Only  carriers  generated  within  a  diffusion  length  of  the  junction 
are  able  to  make  it  to  the  junction  in  time  before  they  are  recombined. 
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Device  Characteristics 


Photoresponse.  A  measure  of  the  performance  of  a  photodetector  is 
its  responsivity.  Responsivity,  or  photoresponse,  is  defined  as  the 
amount  of  photocurrent  generated  per  unit  of  optical  power  incident  on 
the  device,  i.e.  responsivity  is  a  measure  of  flk  of  photocurrent  per 
of  optical  power: 


where 

R  =  responsivity,  nk/  u\J 
1  =  photocurrent,  /.'A 
H  =  optical  power,  ti W 

Both  the  current  and  the  power  are  functions  of  wavelength.  The  develop¬ 
ment  of  this  relatively  simple  equation  can  be  explained  using  Figure  8, 
which  shows  the  energy  band  diagram  of  a  p-n  junction  under  reverse  bias. 
The  photocurrent  generated,  assuming  total  collection  of  the  incident 
light  and  generated  carriers,  is 


I 


(8) 


where 

T]  -  quantum  efficiency,  electrons/photon 

optical  power  incident  on  the  surface,  u W 
q  =  charge  of  an  electron, coul 

The  optical  power  decays  exponentially  into  the  device,  similarly  to 
Equation  (5) : 


H(x)  “  H^exp  ("Ctx) 


(9) 


19 


Fig  8.  The  energy  band  diagram  of  a  p-n  junction  under  reverse  bias. 

Recall  that  there  are  window  losses  and  recombination  losses  incurred 
by  the  photodetector. 

If  collection  takes  place  between  x^  and  then 

I  -  q7?H(Xl)  -  H(x2) 
h  V 


and  by  the  definition  of  responsivity : 


I 

H 


exp(-dx^)  -exp(-ax2)J  fjA/  j/W 


(ID 


The  above  expression  reaffirms  the  notion  of  losses  in  a  photodetector 
and  that  the  photoresponse  is  a  measure  of  the  optically  generated  cur¬ 
rent  in  the  depletion  region  of  a  p-n  junction. 

A  parameter  which  influences  the  photoresponse  measurement  is  the 
reflectivity  of  the  material.  In  the  case  of  no  antireflection  coating 

H'  =  (l-fl)H  (12) 


where 

G(  =  reflection  coefficient 
H'  =  incident  optical  power  after  reflection 
The  reflection  coefficient  is  determined  by  (Ref  53:52): 

<?  -  ^=§4  <»> 

(l+n) 

where  n  is  the  refractive  index  (for  Ge,  n  =  3.46). 

Quantum  Efficiency.  Another  measure  of  how  effective  the  junction 
is  in  converting  incident  photons  to  electrons  is  the  quantum  efficiency, 
denoted  by  Y],  and  expressed  in  units  of  electrons  per  photon.  Like  res¬ 
ponsivity,  quantum  efficiency  is  affected  by  depletion  layer  width.  In 
this  case,  it  will  be  determined  by  the  expression  for  responsivity 

R  =  — ^  A/W  (14) 
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Fig  9.  Ideal  current-voltage  characteristics 

(a)  linear  plot 

(b)  semilog  plot  (Ref  53:101). 

Solving  for  T)  : 

T)  =  Ihc  electrons/photon  (15) 

A  q  H  A  q 

Diode  Characteristics 

The  current-voltage  characteristics  of  a  diode  are  shown  in  Figure  9 
for  an  ideal  diode  (Ref  53:101).  The  dashed  line  signifies  the  effect 
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of  light  incident  on  a  photodiode.  The  generation  of  carriers  in  the 
depletion  region  as  a  result  of  photo-illumination  increases  the  reverse 
leakage  current.  The  number  of  carriers  generated  is  directly  propor¬ 
tional  to  the  varying  reverse  leakage  current.  In  terms  of  current  den¬ 
sity,  the  total  current  density  through  the  depletion  region  under  reverse 
bias,  neglecting  thermal  reverse  current  at  steady  state  is 


J  _  =  J.  +  J, ... 
tot  dr  diff 


(16) 


where 


J 


dr 


Jdiff 


=  drift  current  density 
=  diffusion  current  density 


The  drift  current  density  is  due  to  tne  carriers  generated  inside  the 
depletion  region  and  swept  out  by  the  electric  field.  The  diffusion 
current  density  is  due  to  the  carriers  generated  outside  the  depletion 
layer  in  the  bulk  of  the  semiconductor  and  diffusing  into  the  reverse 
biased  junction. 

The  contribution  from  J,  is  determined  by  the  rate  at  which  elec- 

dr 

tron-hole  pairs  are  generated  by  the  incident  photon  flux.  The  gener¬ 
ation  of  the  carriers  drops  off  exponentially  as  expected  because  light 
intensity  is  absorbed  exponentially  as  a  function  of  distance  into  the 
semiconductor. 

For  distances  outside  the  depletion  layer,  the  minority  carrier 


< 


density  (electrons)  in  the  bulk  of  the  semiconductor  is  obtained  from 
the  one-dimensional  diffusion  equation  (Ref  20:126;  53:99): 


(17) 


2 

D  d  n  n  -n 
n  p  _  p  po 


d  x 


r 

n 


=  0 


where 


2 

=  diffusion  coefficient  for  electrons,  cm  /sec 
=  lifetime  of  excess  carriers,  sec 

_3 

=  density  of  electrons  at  equilibrium,  cm 

_3 

=  minority  carrier  density  (electrons) ,  cm 


Equation  (17)  is  subject  to  the  boundary  condition 


n 

=  n 

for  x  =  ® 

P 

po 

n 

=  0 

for  x  =  W 

P  . 

which  leads  to  the  following  solution  for  total  current: 


(18) 


J  =  Js  [exp(qV/kT)  -l]  (19) 

Under  reverse  bias,  the  current  density  (Ref  53:100)  should  saturate  at 


J 

s 


qD  p 
p  no 


U  n 
q  n  po 


(20) 


where 


L 

p,n 


2 

=  diffusion  coefficient  of  holes,  cm  /sec 

-3 

=  density  of  holes  at  equilibrium,  cm 

■  diffusion  length  of  holes  and  electrons,  respectively, 


H  m 
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For  Ge,  the  intrinsic  carrier  concentration,  n^  =  2.4  x  10  cm 

15  -3  -3 

and  for  3  Q- cm  substrate  wafers,  N  =  1.1  x  10  cm  .  Using  =  10 

seconds  (Ref  53:58)  as  the  minority  carrier  lifetime  and  kT/q  =  0.0259  eV, 

all  values  in  Equation  25  are  known  except  for  the  diffusion  coefficient, 

D  : 
n 


(26) 


2 

where  fi  is  the  drift  mobility  of  electrons,  3900  cm  /V-sec  for  Ge 
(Ref  53:58).  Appropriate  substitution  into  Equation  25  leads  to  an  est¬ 
imate  value  for  the  saturation  current  density: 

Js  =  2.67  x  10"5  A/ cm2  (27) 

To  get  an  idea  of  the  reverse  leakage  current  that  can  be  expected  from 

shallow  junction  photodiodes,  the  current  density  must  be  multiplied 

by  the  area  of  the  junction  active  in  photogenerated  current. 

Dark  current.  When  no  light  is  incident  on  the  photodiode,  the 

current  produced  by  the  diode  is  known  as  the  dark  current  of  the  detector, 

denoted  by  I,.  The  dark  current  is  a  result  of  bulk  and  surface  effects, 
a 

The  bulk  leakage  current  depends  on  the  active  area,  the  resistivity  of 
the  substrate,  and  the  bias  voltage.  The  surface  leakage  current  is  a 
result  of  defects  on  the  surface  of  the  diode,  which  depends  on  cleanli¬ 
ness  during  the  fabrication  of  the  diode.  It  is  also  a  function  of  bias 
voltage  and  surface  area. 

The  dark  current  of  a  photodetector  reduces  the  sensitivity  of  the 
device  because  it  produces  shot  noise  as  given  by 

!  .  =  (2ql  B)**  (28) 

noise  d 

where 


I  .  =  noise  current 

noise 

B  =  noise  equivalent  bandwidth  (Hz) 

The  shot  noise  produced  by  the  dark  current  dominates  the  Johnson  (thermal) 
noise  of  the  diode  in  this  case. 
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Some  methods  of  reducing  the  dark  current  can  be  used  during  the 
design  and  fabrication  process.  A  deeply  diffused  guard  ring  around  the 
active  area  of  the  diode  helps  reduce  the  dark  current  by  shunting  the 
surface  leakage  current  away  from  the  load  resistor.  Again  cleanliness 
during  device  processing  is  a  very  important  factor  in  reducing  the 
leakage. 

Carrier  lifetimes.  Note  that  the  Equation  (25)  expresses  the  sat¬ 
uration  current  density  as  a  function  of  carrier  lifetime,  T  .  A  longer 
lifetime  implies  reduced  leakage  current,  and  greater  sensitivity.  The 
assumption  of  7’n~  1  millisecond  turns  out  to  be  rather  high  when  compared 
to  the  measured  lifetimes  of  the  devices  studied  in  this  thesis  as  will 
be  seen  later.  Equation  (24)  shows  the  direct  relationship  between  dif¬ 
fusion  length  and  longer  lifetime.  A  longer  diffusion  length  allows 
detection  at  longer  wavelengths. 

Applied  Bias.  The  effect  of  applied  bias  was  discussed  earlier  in 
this  section  in  conjunction  with  Figures  5  and  6.  It  can  be  seen  from 
this  figure  and  Equation  (6)  that  the  depletion  layer  width,  and  in  turn 
the  sensitivity,  increases  with  applied  bias  voltage.  Of  course,  at  a 
certain  point,  the  increase  in  response  is  outweighed  by  increased  noise. 

Fabrication  of  Shallow  Junctions 

The  fabrication  of  shallow  junction  Ge  photodiodes  can  be  done  by 
two  well  known  processes,  diffusion  and  ion  implantation.  Although  dif¬ 
fused  diodes  promise  uniform  doping  across  the  surface  of  the  wafer,  the 
depth  of  the  diffusion  is  difficult  to  control  for  shallow  junctions. 
Diffusion  of  impurities  in  a  Ge  wafer  by  an  infinite  source  is  modeled 
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by  the  error  function  complement  (Ref  20:46;  53:82): 


C(x,t)  =  C  erfc  — —  1  cm 

S  2-^Dt  J 


(29) 


where 

C(x,t)  =  impurity  concentration  as  a  function  of  time  and  distance 
from  the  surface,  cm 

-3 

=  impurity  concentration  at  the  surface,  cm 
x  =  distance,  cm 

2 

D  =  diffusion  coefficient,  cm  /sec 
t  =  time,  sec. 

Cg  is  fixed  by  the  solid  solubility  of  the  impurity  in  Ge.  For  a  phos¬ 
phorus  diffusion  at  750°C,  the  solid  solubility  of  phosphorus  in  Ge  is 
20  3 

~ 2x10  atoms/cm  and  the  corresponding  diffusion  coefficient  is 
-12  2 

~  2.8  x  10  cm  /sec  (Ref  53:31,32).  It  should  be  noted,  however,  that 
in  practice  diffusion  coefficients  for  most  impurities  vary  significantly 
from  furnace  to  furnace  probably  because  of  temperature  variations.  The 
impurity  concentration  and  junction  depth  of  shallow  junction  diodes  can 
be  better  controlled  by  ion  implantation.  Implanted  diodes  studied  in 
this  thesis  were  fabricated  using  As  as  the  impurity  species.  Details 
will  follow  in  Chapter  IV. 

In  considering  the  diode  from  the  perspective  of  a  photodetector 
application  edge  breakdown  occurs  more  easily  in  shallow  junctions 
(Ref  20).  There  are  various  methods  of  reducing  edge  breakdown  effects 
which  involve  guard  ring  structures,  diffused  pockets  of  the  window 
material  on  the  outside  of  the  shallow  junction  active  area.  Guard  rings 
encourage  uniform  breakdown  over  the  active  area  of  the  diode,  assuming 
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there  are  no  defects  in  the  material  called  microplasmas,  and  help  reduce 
the  leakage  current  by  shunting  it  to  ground. 

The  surface  leakage  currents  associated  with  shallow- junction  Ge 
planar  diodes  can  also  be  reduced  by  forming  mesa  structures.  These  are 
formed  by  etching  the  surrounding  surface  area  of  a  diode  to  a  depth  below 
its  junction.  However,  mesa  structures,  which  are  an  absolute  necessity 
in  photodiodes  made  from  diffused  wafers,  have  fabrication  problems. 

Based  on  practical  experience,  the  etching  process  is  a  "dirty"  step  and 
is  less  convenient  than  planar  processing.  Consequently,  implanted 
photodiodes  suggest  easier  processing  applications  because  they  have  no 
mesa  etching  requirement. 

Passivating  films  can  also  be  employed  to  reduce  surface  leakage 
currents;  however,  there  is  no  workable  oxide  for  Ge  that  is  analogous 
to  its  Si02  couterpart  that  has  been  found  to  be  so  useful  and  versatile 
in  Si  technology.  The  Ge02  is  water  soluble  and  highly  volatile  so 
Si^N^  is  being  tried.  Various  combinations  of  these  three  methods: 

«uard  rin^s,  mesas,  and  nassivatino  films,  can  be  emnloved  to  reduce  the 
leakage  currents  which  affect  the  photoresponse  of  the  Ge  photodiodes. 

There  are  many  possibilities  for  making  contact  to  the  diode.  One 
method  is  demonstrated  in  Figure  10(a).  A  concentric  ring  of  metalliza¬ 
tion  overlaps  the  active  area  of  the  diode  and  contact  is  made  to  a  "tab" 
extending  from  the  ripg  on  the  outside.  This  technique  requires  an  in¬ 
sulating  oxide  as  shown.  The  simpler  technique  shown  in  Figure  10(b) 
is  used  for  the  diodes  studied  in  this  thesis.  A  concentric  center  con¬ 
tact  does  not  require  the  use  of  an  oxide  but  the  pressure  of  the  contact 
lead  on  the  shallow  junction  can  cause  problems.  It  is  assumed  that  no 
optical  radiation  penetrates  the  metallization,  therefore  the  active 


Fig  10.  Making  contact  to  the  photodiode. 

area  of  the  diode  is  the  total  area  minus  the  contact  area. 

In  summary,  the  operation  of  a  photodetector  is  based  on  the  tran¬ 
sition  of  carriers  from  the  valance  band  to  the  conduction  band  of  the 
semiconductor  material.  The  transition  occurs  by  means  of  absorption 
of  the  photon  if  optical  radiation  whose  energy  must  be  greater  than  or 
equal  to  the  energy  gap  of  the  semiconductor.  Photoexcited  carriers  are 
collected  in  the  depletion  region  of  a  p+n  junction.  The  photocharacter¬ 
istics  of  a  photodiode  are  described  in  terms  of  responslvity,  quantum 
efficiency,  and  dark  current  along  with  known  diode  performance  character¬ 
istics.  The  shallow  junction  structure  is  intended  to  increase  these 
performance  ratings  by  collecting  a  greater  percentage  of  light  absorbed 
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by  the  material.  Specific  fabrication  steps  for  five  shallow- junction 
Ge  photodiodes  are  discussed  in  Chapter  4.  Immediately  following  this 
chapter  is  a  theoretical  discussion  of  the  Ge-Si  alloy  of  the  future 
whose  energy  band  structure  resembles  that  of  pure  Ge.  This  explains 
why  fabrication  proceedures  for  the  photodetectors  of  interest  reflect 
the  material  properties  of  pure  Ge. 
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Ill  Ge-Si  Alloy  Theory 


Two  areas  of  research  are  presently  in  progress  at  the  Deputy  for 
Electronic  Technology  of  the  Rome  Air  Development  Center  relevent  to  the 
goal  of  photodetection  near  1.27  The  area  in  which  research  is  pro¬ 

gressing  quickly  is  that  of  device  fabrication  and  testing  of  shallow- 
junction  structures,  it  is  possible  to  obtain  higher  responsivity  and 
efficiency  than  those  of  comparable  devices  having  deeper  junctions. 

The  second  area  of  research  is  the  crystal  growth  of  the  Ge-Si  alloy 
which  is  intended  to  provide  the  approximate  value  of  the  energy  bandgap 
required  for  optimum  photodetection  near  1.27  u m,  i.e.,  the  region  of 
minimized  attenuation  and  dispersion  in  an  optical  fiber.  Eventually 
these  two  lines  of  work  will  join  in  a  single  effort  to  fabricate  shallow- 
junction  photodetectors  using  the  Ge-Si  alloy  as  the  starting  material. 
Once  a  working  device  has  been  achieved,  it  will  be  tested  to  determine 
if  indeed  the  Ge-Si  photodetector  exhibits  peak  photoresponse  near  1.27  /m 
There  is  a  third  area  of  research  which  should  be  initiated  in  par¬ 
allel  with  the  two  mentioned  above.  The  physics  involved  in  the  binary 
alloy  must  not  be  neglected.  A  quantum-mechanical  study  of  the  resulting 
energy  band  structure  of  the  Ge-Si  system  would  be  desirable  as  an  aid 
in  interpreting  experimental  results.  Theoretical  calculations  of  energy 
bands  have  significantly  predicted  the  qualitative  electrical,  optical, 
and  thermal  properties  of  semiconductors,  although  quantitative  properties 
must  generally  be  obtained  experimentally  (Ref  8) .  Theoretically,  the 
band  structure  of  a  crystal  is  described  by  its  energy -momentum  (€-k) 
relationship. 
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In  this  section,  a  brief  description  of  the  methods  of  determining 
energy  band  structures  will  be  given.  They  include  the  free-electron 
approximation,  the  tight-binding  approximation,  the  orthoganalized  plane 
waves  (OPW)  method,  the  method  of  pseudopotentials,  and  others.  The  LCAO 
method  of  the  tight-binding  approximation  is  chosen  as  a  likely  approach 
for  obtaining  the  energy  band  structure  of  the  Ge-Si  alloy.  An  outline 
of  the  LCAO  is  given  and  applied  to  the  Ge-Si  alloy  in  terms  of  desired 
parameters . 

The  One-Electron  Model 

This  section  begins  with  the  one-electron  model  to  introduce  elemen¬ 
tary  band  theory.  The  one-electron  model  has  been  used  to  develop  the 
theory  of  energy  levels  in  crystals.  The  model  leads  to  a  qualitative 
energy-level  structure  which  is  similar  to  those  in  real  crystals  and 
provides  and  understanding  of  the  behavior  of  electrons  under  the  influ¬ 
ence  of  applied  electric  fields.  The  wave  nature  of  electrons  suggests 
a  quantum-mechanical  analysis  to  obtain  a  proper  description  of  electron 
behavior  (Ref  8:2-3;  15:492-496;  23:1703-1710;  35:208-224;  44:3-11; 
49:196-205). 

The  equation  for  an  electron  of  mass,  m  moving  in  a  region  (the 
crystal  lattice)  whose  potential  energy  is  U(x),  is  the  wave  equation 
known  as  Schroedinger 's  equation: 

-*L.  d2#  +  uW  &  -  Jlk  d#  (30) 

2m  2  d  t 

dx 

This  is  a  partial  differential  equation,  for  the  moment  reduced  to  one 
dimension,  where  the  two  indenendent  variables  are  oosition,  x,  and  time, 
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The  solution  is  the  wave  displacement  or  wavefunction,  denoted  by 
which  is  the  dependent  variable.  If  a  particle  is  in  a  stationary  state, 
its  wavefunction  is  independent  of  time  and  dependent  on  position  only. 

It  can  be  factored  into  separate  position  and  time  factors  and  is  expres¬ 
sed  by 

#(x,t)  =  0(x)e"jfc/'fi)t  (31) 

where  €.  is  the  total  energy  of  the  electron,  potential  plus  kinetic. 
Substitution  of  this  proposed  solution  into  the  wave  equation  produces 
another  differential  equation  for  ^(x)  known  as  Schroedinger 's  time 
independent  equation: 

-*2  ^  +  u(x)^  =€</>.  (32) 

2m  >2 
a  x 

Two  factors  have  been  neglected  which  greatly  complicate  the  quan¬ 
tum  mechanical  problem.  First,  each  electron  is  affected  by  a  combination 
of  forces  produced  by  all  the  other  electrons.  Second,  the  thermal  vi¬ 
bration  of  the  lattice  occurs  in  many  different  modes,  and  the  resulting 
time-dependent  forces  are  also  applied  to  all  the  electrons  thru  the 
potential  term,  U(x).  The  one-electron  model  assumes  that  the  motion  of 
an  electron  is  independent  of  all  the  other  electroris  and  of  the  lattice 
vibrations  so  that  the  environment  of  the  electron  is  that  of  a  strictly 
periodic  potential.  Under  these  assumptions,  the  solution  of  the 
Schroedinger  equation  for  the  single  electron  is  a  set  of  states  which 
can  be  occupied  by  the  single  electron  or  any  other  electron  in  the 
crystal.  The  solutions  are  one-electron  states,  however,  limited  by  the 

( 
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Pauli  exclusion  principle  which  says  that  only  one  electron  can  occupy 
a  given  energy  state  at  one  time. 

A  periodic  crystal  potential  in  one  dimension  along  one  of  the  (100) 
directions  in  a  cubic  crystal  is  shown  in  Figure  11  where  a  is  the  lattice 
constant.  The  high  surface  potential  of  a  crystal  is  neglected  in  favor 
of  extending  the  periodicity  of  the  crystal  infinitely  in  all  directions, 
or  forming  a  closed  ring  with  the  periodic  potential.  The  solutions  to 
the  one-electron  model  are  wavefunctions  called  Bloch  functions 

^(x)  =  e*^  u(x)  (33) 

The  properties  of  these  wavefunctions  are  closely  related  to  the  period¬ 
icity  of  the  lattice,  and  the  allowed  energy  states  of  the  electrons 
form  continuous  bands  of  allowed  states  which  are  separated  by  forbidden 
energy  regions. 

The  Kronig-Penney  model  of  an  infinite  one-dimensional  crystal  as¬ 
sumes  the  infinite  one-dimensional  square-well  potential  shown  in 
Figure  12.  Although  it  is  a  crude  approximation  compared  to  an  actual 
crystal,  it  is  useful  in  demonstrating  some  of  the  important  character¬ 
istics  of  the  quantum  behavior  of  electrons  in  a  periodic  lattice.  To 
solve  the  one-electron  model  using  the  square-well  potential, 

Schroedinger 's  equation  takes  the  form 

%  +  -2?  [€  -V(x)]  ^(x)  =  0.  (34) 

a  x  h 

The  wavefunctions  must  have  the  Bloch  form  in  each  interval  of  the  per¬ 
iodic  potential,  i.e.  where  the  potential  is  high  and  where  it  is  low. 
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Fig  11.  A  representation  of  the  potential  in  a  perfectly  periodic 

crystal  lattice.  The  surface  potential  barrier  is  at  the  left 
(Ref  35:209). 


Fig  12.  The  ideal  periodic  square-well  potential  used  by  Kronig  and 

Penny  to  show  the  general  characteristics  of  the  quantum  behavior 
of  electrons  in  periodic  lattices  (Ref  35:213). 
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Fig  13.  The  energy  €  plotted  as  a  function  of  k  (Ref  35:217). 

The  resulting  differential  equations  obtained  by  substituting  the  proposed 
solutions  into  the  Schroedinger  equation  are  solved  by  standard  procedures 
and  appropriate  boundary  conditions  are  applied.  Legitimate  solutions 
are  those  which  produce  well-behaved  wavefunctions  which  are  physically 
realizable.  They  are  those  which  satisfy  the  physical  requirements  of 
the  system. 

A  plot  of  energy  €  as  a  function  of  k  as  shown  in  Figure  13  reveals 
alternating  regions  of  allowed  and  forbidden  energy  bands  characteristic 
of  the  behavior  of  electrons  in  periodic  lattices.  For  large  energies, 


€  (k)  approaches  the  free  electron  relation  (Ref  35:79-80) 


6 


2  2 

=  fivr 

2m 
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shown  by  the  dotted  line  within  the  allowed  bands.  Here,  the  allowed 
bands  become  very  broad  compared  to  the  narrower  forbidden  bands. 

When  an  electric  field  is  applied  to  the  crystal,  the  momentum  of 
the  electron  in  the  periodic  lattice  changes.  The  effect  of  the  periodic 
crystal  potential  is  to  replace  the  mass  of  a  free  electron  with  an 
effective  mass  described  by 


m* 


(36) 


which  plays  the  role  of  the  mass  in  Newton's  force  equation.  When  the 
€  -k  curve  is  parabolic,  which  the  curve  in  Figure  13  approximates,  the 
effective  mass  becomes  a  constant.  This  is  a  good  first  approximation 
of  the  effect  of  the  periodic  crystal  potential  on  a  free  electron.  The 
consequence  is  that  in  equations  involving  the  mass  of  an  electron,  the 
free-electron  mass  must  be  replaced  by  an  appropriate  effective  mass 
which  is  related  to  the  form  of  the  €-k  curve  as  in  Equation  (36) .  For 
example,  the  mobility  fl  of  a  carrier  is  a  function  of  the  effective  mass 
(Ref  22:165)  by 


5JL 

2  m* 


(37) 


and  is  thus  related  to  the  energy  band  structure  through  the  slope 
3  €  / 5  k.  It  should  be  made  clear  that  the  wavefunctions  found  as 
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Fig  14.  The  reduced-zone  representation  of  the  f -k  plot  of  Figure  13 
(Ref  35:220). 

solutions  to  the  wave  equation  actually  apply  to  a  wavepacket  of  many 

electrons  rather  than  one  electron.  Therefore  the  effective  mass  applies 

to  a  wavepacket  rather  than  one  electron.  Fig  13  can  be  conveniently 

condensed  as  in  Fig  14  by  translating  all  segments  of  the  versus  k 

curve  into  the  first  Brillouin  zone,  -  —  <  k  <  where  a  and  b 

a+b  ~  ~  a+b 

are  defined  by  Figure  12. 

The  one-electron  model  with  a  periodic  potential  has  so  far  been 
used  to  obtain  the  energy  band  structure  of  a  crystal  (€  -k) ,  demonstra¬ 
ting  its  similarity  to  the  free-electron  model  in  terms  of  allowed  and 
forbidden  energy  states.  It  also  leads  to  the  modification  of  the  free- 

electron  model  by  the  replacement  of  the  free  electron  mass  by  its 
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effective  mass,  a  funtion  of  the  slope  of  the  €  -k  curve  in  the  allowed 
energy  bands.  The  model  was  based  on  the  assumption  that  the  motion  of 
an  electron  was  independent  of  all  the  other  electrons  and  of  the  lattice 
vibrations,  dependent  only  on  the  periodic  potential  of  the  crystal.  In 
actual  crystals,  the  potential  function  which  influences  the  motion  of 
and  electron  must  be  related  to  the  ion  cores  and  the  other  electrons  of 
the  crystal.  An  exact  solution  to  this  problem  is  impossible;  therefore, 
the  problem  must  be  approached  using  an  appropriate  approximation. 

The  Free-Electron  Approximation 

The  free-electron  approximation  is  based  on  the  assumption  that  the 
total  energy  of  the  electron  is  always  large  compared  to  the  energy  of 
the  periodic  potential.  The  resulting  allowed  bands  will  be  broad  and 
the  forbidden  energy  regions  will  be  very  narrow  (Ref  35:224-231;  44:10-15 
49:205-208).  In  actual  crystals,  the  periodic  potential  goes  to  -oo  at 
the  ion  cores,  but  this  approximation  holds  fairly  well  for  the  outermost 
electrons  in  metals  and  alkali  metals. 

The  first  step  is  to  represent  the  periodic  potential  V(x)  in  a 
simplified  one-dimensional  model  as  (Ref  35:225): 

=  yf(x)  (38) 
n 

where 

m  =  free  mass  of  an  electron 
h  =  Planck’s  constant  divided  by  27T,  h/27T 
y  **  scalar 

and  f(x)  has  the  periodicity  of  the  lattice.  The  Schroedinger  equation 
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takes  the  form 


+ 


k  + 
o 


where  kQ  is  related  to  the  total  energy  as  follows 


2  2 
ft  k 

o 

2m 


(39) 


(40) 


In  solving  the  Schroedinger  equation,  the  periodicity  of  f(x)  allows  it 
to  be  expressed  as  a  Fourier  series  whose  coefficients,  of  f(x),  are 
also  related  to  the  Fourier  expansion  coefficients  of  the  periodic 
potential  V(x).  The  wavefunctions  take  the  Bloch  form  and  are  in  part, 
through  u^(x) ,  a  periodic  function  of  period  a,  and  can  thus  also  be 
written  as  a  Fourier  series. 

In  the  case  of  a  free  electron,  y  =  0  and  the  solutions  are  the 
wavefunctions  of  a  free  particle  (Ref  35:226): 

ik  x 

<^(x)  =  bQe  °  (41) 

where  kQ  is  defined  in  Equation  (39)  and  bQ  is  a  constant.  This  solution 
can  also  be  written  as  a  Fourier  series  such  as  (Ref  35:226): 

<//(x)  -  b  eikx  +  y  Likx  Sbe-27r  in  x/al 

°  [  n^O  n  J  (42) 


where  n  «  1,2,3,.., 

k  **  propagation  constant 
a  «■  lattice  constant 

and  whose  coefficients  approach  zero  except  for  bQ.  For  small  values 
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of  y  specified  by  the  f ree-electron  approximation  the  solution  is  valid, 
The  solution  is  recognized  to  have  the  Bloch  form  where 

u  (x)  -  b  +  vL  e2  TTinx/a 
k'  o  n*0  n  (43) 


I  XcCX 

is  obtained  by  dividing  t//(x)  by  e  .  The  solution  is  substituted  into 
the  Schroedinger  equation  producing  a  quadratic  equation  in  y  set  equal 
to  zero.  After  lengthy  manipulation  (Ref  35:226-229),  the  6-k  relation¬ 
ship  becomes  (Ref  35:229): 

^  ♦  (*-  +(^N)2  (44) 


At  the  band  edges  k  =  k^  =  n7T/a,  6  (k)  becomes 


€  =  £  +  V 

n  -  n 


where 


f„  ■  <« 

represents  the  energy  of  a  free  particle  associated  with  the  band  edge 


points  and 


V  =^c 
n  2m  n 


where  is  the  nth  Fourier  coefficient  in  the  expansion  of  the  periodic 
potential  as  a  Fourier  series  (Ref  35:229). 

At  the  band  edges,  the  discontinuity  is  due  to  internal  Bragg  reflec¬ 
tion.  The  £-k  curves  shown  in  Figure  15  show  that  the  curves  are  nearly 
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( 


Fig  15.  Schematic  representation  of  the  €  versus  k  relation  for  the 
free  electron  approximation  (Ref  35:230). 


parabolic  at  the  band  edge  points.  This  complicates  the  differentiation 
necessary  to  determine  the  effective  mass,  so  the  calculation  is  simpli¬ 
fied  by  using  another  approximation 


k 


n  IT 
a 


+  k' 


where  k’ «  7T /a  which  is  substituted  into  Equation  (44).  The  width  of 
the  forbidden  energy  region  at  the  band  edge  k  =  n  7T/a  is 

A€  =  2  |vn  |  (48) 

and  with  the  aid  of  the  binomial  theorem  (Ref  35:230-231): 

,2.  ,2 

€00  -  €n  +  *iA€  +S2S~  (49) 
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The  effective  mass  becomes  (.Ref  35:231): 


where 

m  =  free  mas  of  an  electron 
£  =  energy  of  a  free  particle  at  the  band  edge 

=  width  of  a  forbidden  energy  region 
which  is  the  same  for  holes  at  the  top  of  the  lower  band. 

The  result  of  the  free-electron  approximation  depicts  an  £-k  curve 
with  broad  allowed  energy  states  and  narrow  forbidden  regions.  Raimes 
develops  the  same  model  by  perturbation  theory  and  use  of  the  Fourier 
expansion  and  explains  the  Bragg  reflection  of  electrons  at  the  band 
edges  (Ref  49:205-208).  Another  way  to  look  at  the  assumption  on  which 
the  free-electron  approximation  is  developed  is  to  say  that  the  atoms  of 
the  crystal  are  so  close  together  that  the  wave  functions  for  electrons 
of  neighboring  atoms  overlap  somewhat.  The  strong  interaction  between 
atoms  implies  that  the  wavefunctions  representing  the  allowed  energy 
states  in  the  lattice  are  not  very  similar  to  the  wavefunctions  describing 
the  individual  atoms  of  the  crystal.  There  is  another  approximation  which 
takes  the  opposite  point  of  view,  discussed  in  the  following  section. 

The  Tight-Binding  Approximation 

The  tight-binding  approximation  takes  the  opposite  point  of  view  to 
the  free-electron  approximation.  In  this  case,  the  approximation  is  based 
on  the  assumption  that  atoms  of  the  crystal  are  far  apart  so  the  inter¬ 
action  between  atoms  is  very  weak.  The  resulting  wavefunctions  and 
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allowed  energy  states  of  the  electrons  will  be  closely  related  to  the 
wavefunctions  and  allowed  energy  levels  of  the  individual  atoms.  The 
results  predict  narrow  allowed  energy  bands  in  comparison  with  the  for¬ 
bidden  bands  (Ref  35:231-236;  49:208-214).  The  particular  crystal  deter¬ 
mines  the  appropriate  approximation  but  it  may  also  be  necessary  to  find 
an  intermediate  solution  for  crystals  which  do  not  fit  in  either  category. 
In  many  cases,  the  tight-binding  approximation  is  best  suited  for  semi¬ 
conductors  and  the  free-electron  approximation  is  more  suitable  for 
metals  (Ref  44:16). 

The  wavefunctions  of  the  crystal  in  the  tight-binding  approximation 
are  based  on  the  wavefunctions  of  the  individual  isolated  atoms.  If 
V  (  F)  represents  the  potential  function  of  an  isolated  atom,  the  solu¬ 
tions  of  the  Schroedinger  equation  will  represent  the  wavefunctions  of 
that  atom.  The  wave  equation  takes  the  form  (Ref  35:232): 

zlL.  F2  +  V  (r)  0(  r)  =  €  0  (7)  (51) 

2m  °  °  vo 

where  4,Q  represents  the  nondegenerate  wavefunction  of  the  ground  state 
corresponding  to  the  energy  of  the  ground  state.  When  the  wave- 

functions  of  the  nearby  lightly  bound  atoms  of  the  crystal  are  brought 
into  the  picture,  the  crystal  wavefunction  becomes  a  linear  superposition 
of  atomic  wavefunctions  (Ref  35:232): 

^(7)  =  ?an  ^o(  7"7n}  *  (52) 

The  vectors  r  and  rR  are  shown  in  Figure  16  and  the  sum  is  taken  over  all 
the  atoms  of  the  crystal,  n.  Since  all  the  atoms  of  the  crystal  are  the 
(  same,  the  coefficients  a^  must  have  the  same  absolute  value  and  must  be 
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Fig  16.  The  r  and  r  vectors  used  in  the  calculation  of  the  tight- 
binding  approximation  (Ref  35:232). 

±d>  , 

able  to  be  expressed  in  the  form  ae  n  where  a  and  <p ^  are  real  values. 
The  phase  difference  between  neighboring  atoms,  denoted  by  </>,  must  be 
the  same  between  all  atoms  of  the  crystal  along  a  specified  crystal  axis. 
A  phase  factor 


< i>  =  k  •  r 

’n  n 


(53) 


is  chosen  such  that  k  is  a  constant  vector,  and  a  is  set  to  1  to  preserve 
the  normalization  of  (Ref  35:232).  The  wavefunctions  take  the  form 


0_(  r) 
k 


'r)  ■  JLj  e*k  *  rn  ^  (r-r  ) 


(54) 
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which  are  in  Bloch  form  and  satisfy  the  Schroedinger  equation  for  the 
entire  crystal.  The  periodic  potential  function  of  Figure  17  is  used  in 
this  equation  as 


zL.  ^ 2  (/>_  +  v(r)  \p_  fjj 

2m  '  V  *k  k  €  ^k 


(55) 


The  Hamiltonian  of  this  equation  is 


"‘Ti  V2  +Vo(r) 


(56) 


and  is  broken  up  into  two  parts  (Ref  35:233): 


H  =  H  +  H' 
o 


(57) 


where 


h2  9 

H  =  =~  V  *  +  V  (r-r  ) 
o  Zm  on 


(58) 


and 


H'  -  V(r) — V  (r-r  ). 

o  n 


(59) 


For  the  ground  state  wavefunction 


H  lb  =  £  lb 
0*0  o  *  o 


(60) 


and 


H0  “  Zyk*rn  H0  ^Q(r-rn) 

h 

"  C„]£eik#rn  ^0(r-rn) 


me  lb 
0  k 


(61) 
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Fig  17.  The  potential  function  used  in  the  tight-binding  approximation 
(Ref  35:233). 
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The  discrete  levels  of  atomic  energy  split  into  bands  as  the 
isolated  atoms  are  assembled  into  a  crystal  lattice  (Ref  35.235) 


which  is  simply  the  Hamiltonian  operator  passed  through  the  linear  super¬ 
position  of  atomic  wavefunctions  in  the  equation  for  lp_( r) .  The  constant 

k 

£  is  passed  through  the  summation  sign  leaving  the  results  shown.  The 
energy  €  can  be  obtained  by  finding  the  expectation  value  of  the  Hamil¬ 
tonian  operator  and  requiring  the  following  conditions: 

1.  Neglect  the  samll  overlapping  that  occurs  between  atomic  wave- 
functions. 

2.  This  implies  the  normalization  properties  of  and 

* 

dv  =  N  (62) 

where  N  is  the  number  of  atoms  in  the  crystal. 

The  derivation  (Ref  35:233-235)  involves  the  volume  integral  of  sums,  in 
which  case  the  property  of  linearity  is  taken  advantage  of.  Decreasing 
contributions  of  the  terms  in  the  final  sum  allow  neglecting  the  contri¬ 
butions  of  terms  from  atoms  which  are  not  nearest  neighbors  of  the  atom 
under  consideration.  The  spherical  symmetry  of  the  wavefunctions  is 
assumed  as  a  result  of  this  approximation.  So  far,  the  case  under  con¬ 
sideration  is  restricted  to  the  ground  state  electronic  configuration  in 
an  isolated  atom,  its  s-state. 

If  the  vector  r  connects  an  atom  at  the  origin  to  its  nearest 
neighboring  atom,  at  the  origin  (Ref  35:234): 

-  Ct  =  ^Q*(r)  V(r)-VQ(r)  ^(r)dv  (63) 

while  for  the  nearest  neighbors  (Ref : 35: 234) : 

-p  -  /  0o*(r-rm)  V(r)-VQ(r)  ^o(r)dv.  (64) 
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The  £  -k  relationship  obtained  by  summing  over  the  nearest  neighboring 
atoms  is  (Ref  35:234): 

€(k)  =  £  -a~B  De"lk,rm  (65) 

o  ^  in 

where  Of  and/? are  determined  from  Equations  (62)  and  (63). 

The  range  of  allowed  energy  values  corresponds  to  the  vector  and 
forms  an  energy  band  in  the  same  way  as  in  the  free-electron  approxima¬ 
tion.  The  minimum  value  of  the  energy  occurs  where  k  is  the  zero  vector 
and  the  maximum  value  of  energy  depends  on  the  lattice  constant  a  which 
appears  in  the  components  of  r  .  In  the  limit  as  the  atomic  spacing 
becomes  large  the  tight-binding  model  becomes  a  very  good  approximation. 
The  allowed  energy  bands  become  narrow  and  approach  a  single  discrete 
energy  value,  the  s-state  of  an  isolated  atom,  as  shown  in  Figure  18. 

An  electron  with  low  momentum  will  he  found  at  the  bottom  of  an 
energy  band  where  £  versus  k  is  a  parabolic  relationship.  The  electron 
then  behaves  like  a  free  electron.  The  effective  mass  of  the  electron 
in  a  particular  direction  is  given  by  the  following  equation  using  the 
slope  of  the  £-k  curve  in  one  dimension  (Ref  35:236): 


where  a  is  the  lattice  constant  and  is  defined  by  Equation  (64).  The 
effective  mass  is  an  isotropic,  that  is,  not  the  same  in  all  directions 
of  the  crystal.  In  the  limit  of  the  tight-binding  approximation,  the 
electron  becomes  heavy  and  not  transferred  easily  from  one  atom  to 
another. 

The  tight-binding  approximation  is  also  known  as  the  LCAO  (Linear 
Combinations  of  Atomic  Orbitals)  approximation  (Ref  52)  and  is  also 


described  in  connection  with  the  hydrogen  molecule  (Ref  49: Section  5.8). 
The  LCAO  method  has  many  qualitiative  features: 

1.  It  gives  solutions  showing  the  correct  symmetry  properties  of 
the  energy  bands. 

2.  It  provides  easy  access  to  solutions  for  energy  bands  of  the 
Brillouin  zone  whereas  other  approximate  methods  become  to 
difficult  at  other  than  certain  symmetry  points  in  the  Brillouin 
zone. 

3.  By  retaining  only  the  terms  which  make  the  method  qualitatively 
correct,  the  LCAO  method  is  a  convenient  interpolation  tool. 

The  Hamiltonian  operator  can  be  represented  in  the  form  of  a  matrix  whose 
elements  are  related  to  the  energy  components  of  the  Bloch  functions  and 
whose  integrals  are  replaced  by  simplifying  disposable  constants.  The 
LCAO  method  for  the  crystal  lattice  having  a  diamond  structure  as  is 
found  in  Ge  is  outlined  by  Slater  and  Koster  (Ref  52:1516-1524)  and  will 
be  discussed  in  detail  later. 

Orthogonalized  Plane  Waves  (OPW) 

In  the  orthogonalized  plane  wave  method  (Ref  44:42-43)  the  Bloch 
function  is  expanded  into  many  plane  waves  to  represent  its  rapid  varia¬ 
tions  near  the  ion  core.  The  orthogonalized  funciton  is  a  linear  comb¬ 
ination  of  plane  waves  and  atomic  wavefunctions  of  the  occupied  states  of 
the  ion  cores.  The  new  states  derived  from  the  OPW  methods  show  that  the 
band  energies  are  very  close  to  the  free  electron  energies  except  for  the 
points  in  k-space  very  near  to  the  boundary  of  a  Brillouin  zone  (Ref  29: 
285). 
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The  ion  cores  are  given  by  atomic  orbital  states  assumed  to  be  ortho¬ 
normal  and  related  to  the  periodicity  of  the  lattice.  The  core  functions 
are  given  by(Ref  44:42): 

tr)  -  _1_  y.  u  i-I  )  (67) 

y?  E» 


which  are  related  to  the  periodicity  of  the  lattice  by 


0_  (r+R  ) 


1  y^e^^n  u.Cr+R.-R  ) 
v/n  Rn  3 


(68) 


Linearity  under  summation  allows  this  to  be  written  as 


A.  .  ik«R.  \  ^  ik«(R  -R, )  —  —  . 

<p  (r+R..)  =  e _ 1  /  ,  e  n  1  u.(r+R.-R  ) 

yF  Rn  3 


k.j 


ik»R,  <f>  . 

■e  1 


(69) 


for  R^  replaced  by  R^-R^  (Ref  44:42).  The  resulting  Bloch  function  must 
have  the  periodicity  of  the  lattice. 

ik#r 

The  next  step  is  to  expand  the  wavefunctions  into  plane  waves,  e  , 
orthogonal  to  the  corresponding  wave  vectors  of  the  wavefunction  according 
to  (Ref  44:43): 


K  ■ 6 


ik»r 


(Nfi 


_Z^-  4>~  (7) 

j  M  k,jv  ; 


(70) 


where  the  orthogonality  is  defined  by 


(r)  ( r)d3r 

k 


0 


(71) 
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When  the  equations  for  <f>_  (r+R  )  and  0  (r)  are  substituted  into  the 

k,J 

equation  for  orthogonality,  two  sums  of  integrals  are  obtained  which  re¬ 
duce  to 


(r-R')u.(r-R  )d 
n  3  n 


3 

r  = 


(72) 


assuming  the  atomic  wave  functions  are  orthogonal  and  do  not  overlap 
(Ref  44:43).  The  coefficients  of  expansion,  become  (Ref  44:43): 


V- 

k,l 


f  ik.r  *  ,3 

Jail  6  ul(r)d  r 

space 


(73) 


The  orthogonalized  plane  waves  can  be  obtained  by  substituting  these  co¬ 
efficients  into  the  equation  above.  When  the  electron  is  far  from  the 
ion  core,  the  OPW  emulate  plane  waves.  Near  the  ion  core,  the  OPW  behave 
like  atomic  orbitals  (Ref  44:43).  The  precision  of  the  OPW  method  depends 
upon  the  accuracy  of  the  assumed  potential,  the  accuracy  of  the  eigen¬ 
states  assumed  for  the  core,  and  the  number  of  terms  admitted  to  the 
Fourier  representations  of  the  crystal  wave  functions.  The  OPW  method 
was  applied  to  Si  and  Ge  by  Herman  as  shown  in  Figure  19.  The  ortho- 
gonalization  of  plane  waves  to  ion  cores  is  difficult  so  a  related  method 
is  discussed  below. 


The  Pseudopotential  Method 

The  pseudopotential  method  shows  the  orthogonal  relationship  between 
plane  waves  and  ion  core  wavefunctions  more  clearly  (Ref  44:44-45)  by 
examining  the  orthogonality  between  the  wavefunction  of  an  electron  in 
the  valence  band  and  the  core  function.  This  relationship  is  expressed  by 


\p_  *=  V_  +  'y  a^  0_ 

ki  ki  kj  kj 
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Fig  19.  The  energy  band  structures  of  a)  Si  and  b)  Ge  (Ref  23:1716). 


where 


i 

a_ 

kj 


\  ki  kj  / 

f  ,*  -3 

-/  <p_  V_  (r)d  r 

7  ki  kj 


(75) 


The  problem  is  in  determining  V _  .  By  determining  the  eigenfunctions 

ki 

for  the  crystal  Hamiltonian  (obtained  by  the  tight-binding  approximation) 
which  represent  the  core  funtion  by 


H  <f>_  =  <p_ 

kj  kj  kj 


(76) 


and 


H  ip_  =  lp_ 

ki  ki  ki 


(77) 


the  equation  becomes 


H  V  _ 
ki 


ak1  (€_  -C_)<t>_ 
3  kj  ki  kj 


€_  r_ 
ki  ki 


The  pseudopotential  is  defined  by 


(78) 


V 

P 


=  y2a-  (e_  -  c_  )^ki 

j  kj  kj  ki  T¬ 


(79) 


or 


(H+V  )  V_  =  e_  V_  . 
P  ki  ki  ki 


(80) 


The  pseudopc  tential  represents  a  repulsive  force  which  distorts  the  act¬ 
ual  crystal  potential.  The  results  of  this  method  for  Ge  and  Si  are  shown 
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Fig  20.  The  energy  band  structures  for  a)  Ge  and  b)  Si  using  the 
method  of  pseudopotentials  (Ref  10:739). 

in  Figure  10  (Ref  10:793;  29:316-317). 

The  Cellular  Method 

The  cellular  method  is  an  old,  specialized  method  for  finding  € (k) 
at  a  point  of  high  symmetry  (Ref  23;  44:45-75).  The  method  involves 
solving  the  Schroedinger  equation  in  a  single  atomic  cell  whose  wave- 
functions  must  satisfy  the  boundary  conditions  in  terms  of  the  derivatives 
given  by  Bloch's  Theorem.  Since  boundary  conditions  are  too  complicated 
in  complex  crystals  and  applicable  only  to  simple  ones,  it  is  assumed 
the  cellular  method  would  not  be  suitable  for  the  Ge-Si  alloy.  Also, 
assumptions  based  on  spherical  symmetry  of  the  lattice  potential  may  not 
be  accurate  in  the  binary  alloy  case. 

Effective  Mass  Approach 

An  adaptation  of  the  cellular  method  to  the  states  at  the  bottom  of 


I 


I 

the  lowest  valence  band  form  the  basis  of  this  approach  (Ref  44:46-47). 

i 

|  Since  the  calculations  are  based  on  the  assumption  of  spherical  symmetry 

I 

to  determine  the  boundary  conditions,  it  is  surmised  that  this  approach 
would  not  be  appropriate  to  the  case  of  the  Ge-Si  alloy  for  the  same 
J  reasons  given  for  the  cellular  method. 

Augmented  Plane  Waves 

In  the  augmented  plane  wave  method,  the  wavefunction  of  an  electron 
is  expanded  in  terms  of  a  combination  of  outer  plane  waves  and  inner 
spherical  waves  (Ref  44:47).  Spherical  assumptions  have  been  previously 
mentioned  as  inadequate  for  the  binary  alloy  so  this  method  is  not  advi¬ 
sable  for  theoretical  determination  of  the  6 -k  structure  of  Ge-Si. 

Ge-Si  Alloy  Considerations 

I  Although  an  alloy  is  not  a  perfect  crystal,  one  can  still  speak  of 

its  energy  band  structure.  The  allowed  bands  are  energy  intervals  in 
which  there  is  a  high  density  of  allowed  states,  and  the  forbidden  bands 
are  energy  intervals  in  which  there  is  a  low  density  of  allowed  states 
(Ref  23:1721-1722).  The  Ge-Si  alloy  is  treated  as  having  an  energy  band 
structure  which  is  very  similar  to  perfect  crystals  and  the  optical  gap 
is  considered  to  be  a  measure  of  the  width  of  a  forbidden  band. 

The  Ge-Si  energy  band  structure  has  been  described  previously  in 
Figure  3  (Ref  27;  23:1722;  8:76).  As  Si  is  added  to  Ge  in  the  0-15% 
range,  the  lowest  conduction  band  moves  away  from  the  highest  valence 
band.  The  position  of  the  minima  along  the  (111)  and  (100)  directions 
also  changes  in  terms  of  the  distance  the  minima  move  away  from  the  edge 
(  of  the  valence  band.  Since  the  (111)  minima  lie  below  the  (100)  minima 
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in  the  range  of  interest,  and  the  (111)  minima  in  the  conduction  band 
moves  away  from  the  valence  band  at  a  faster  rate  than  the  (100)  minima, 
the  optical  gap  is  determined  by  the  position  of  the  (111)  minima.  The 
conduction  band  minimum  at  the  origin  moves  away  from  the  valence  band 
edge  much  more  rapidly  than  the  (111)  and  (100)  minima  and  therefore  does 
not  influence  the  measurement  of  the  energy  gap.  Magnetoresistance 
measurements  have  been  shown  to  be  consistent  with  the  theory  thus  des¬ 
cribed  (Ref  23:1722)  for  the  range  of  interest  of  MOL  %  Si  in  Ge. 

The  Tight-Binding  (LCAO)  Approach  to  the  Ge-Si  Alloy 

The  decision  to  apply  the  tight-binding  approximation  to  the  Ge-Si 
alloy  was  based  on  the  approximation's  suitability  for  semiconductors. 

The  equivalent  simplified  LCAO  method  for  periodic  potentials  has  been 
described  by  Slater  and  Koster  for  diamond  structures,  specifically  Ge, 
and  compared  with  Herman's  OPW  calculations  (Ref  52:1516-1524). 

General  Approach  of  the  LCAO  Method  for  Solids.  The  following 
discussion  is  taken  after  a  rather  complete  article  by  Slater  and  Koster 
describing  the  application  of  the  LCAO  method  for  the  periodic  potential 
problem  (Ref  52) .  Recall  that  the  LCAO  method  used  to  solve  the  problem 
of  the  periodic  potential  consisted  of  taking  a  linear  combination  of 
atomic  orbitals  of  the  various  atoms  of  the  crystal.  The  coefficients 
of  the  linear  combination  are  equivalent  to  the  values  of  the  plane  wave 
e  at  the  various  positions  R  where  the  atoms  are  located.  The  method 
can  involve  many  complicated  integrals  which  must  be  evaluated.  These 
integrals  are  treated  as  disposable  constants  which  in  turn  are  chosen 
to  fit  the  results  of  more  accurate  calculations,  such  as  the  OPW  method. 
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at  restricted  points  of  symmetry  in  the  Brillouin  zone  (Kef  52:1498-1502). 


Consider  an  atomic  orbital  </>  (r-R. )  located  on  an  atom  at  a  vector 

n  1 

position  R^,  where  the  subscript  n  denotes  the  quantum  numbers.  The 
Bloch  sum  for  one  atom  in  the  unit  cell  of  the  crystal  (Ref  52:1498): 


E  (v 

position 
of  an  atom 


expCik^R^) 


the  plane  wave 
at  R. 


*„tr-y 

nth  atomic  orbital 
of  an  atom  at  R^ 


(81) 


will  be  extended  over  the  atoms  in  equivalent  positions  in  all  the  unit 
cells  of  the  crystal.  Normalization  will  be  set  up  later.  Such  a  Bloch 
sum  can  be  written  corresponding  to  each  atomic  orbital  of  a  particular 
atom,  and  for  each  atom  in  the  unit  cell  of  the  crystal.  The  set  of  lin¬ 
ear  combinations  is  equivalent  to  the  set  of  sums  of  atomic  orbitals 
operated  on  by  the  Hamiltonian  operator.  The  matrix  of  the  Hamiltonian 
operator  has  no  diagonal  components  representing  the  kinetic  energy  and 
periodic  potential  functions  for  Bloch  sums  with  different  k's.  For 
Bloch  sums  with  the  same  k  values,  there  are  matrix  components  of  the 
Hamiltonian  which  do  not  vanish.  These  correspond  either  to  different 
atomic  orbitals  on  the  same  atom,  or  to  atomic  orbitals  on  different 
atoms  in  the  unit  cell.  Some  nondiagonal  components  of  the  Hamiltonian 
matrix  vanish  at  special  symmetry  points  of  the  Brillouin  zone  or  special 
values  of  k.  The  vanishing  feature  of  some  of  the  matrix  elements  deter¬ 
mines  the  special  properties  of  the  energy  bands  at  the  symmetry  points. 

An  approximate  solution  of  the  periodic  potential  problem  can  be 
set  up  using  the  following  steps  (Ref  52:1498-1500): 

1.  A  finite  set  of  atomic  orbitals  is  taken  on  each  of  the  atoms 
of  the  unit  cell,  starting  with  the  lowest  atomic  states  and 


ending  with  the  highest  levels  occupied  by  the  crystal  Cor  a 
little  higher) . 

2.  A  Bloch  sum  is  constructed  from  each  of  the  atomic  orbitals. 

3.  For  a  given  value  of  k,  a  wave  function  is  set  up  as  a  linear 
combination  of  all  the  Bloch  sums. 

4.  There  are  matrix  components  of  energy  between  all  these  Bloch 
sums  producing  a  secular  problem  of  an  order  equal  to  that  of 
the  total  number  of  atomic  orbitals  considered.  At  special 
values  of  k,  the  secular  equation  can  be  factored  and  thus 
greatly  simplified. 

5.  The  secular  equations  can  be  solved  using  modern  digital  com¬ 
puters  to  give  wavefunctions  and  energy  levels  which  form  use¬ 
ful  approximations  to  the  periodic  potential  problem. 

6.  For  a  given  value  of  k,  the  secular  equation  has  distinct  roots. 
Wavefunctions  of  different  k  values  which  have  the  lowest  energy 
are  grouped  together  into  the  lowest  energy  band,  and  similarly 
for  the  second  energy  band  on  up. 

7.  At  special  symmetry  points,  the  energy  bands  can  be  degenerate. 
The  LCAO  method  described  above  is  a  good  approximation  to  the  periodic 
potential  problem  but  it  requires  a  great  amount  of  numerical  work  to 
compute  the  various  matrix  components  of  energy.  These  will  be  simpli¬ 
fied  later  on.  The  choice  of  one  set  of  atomic  orbitals  over  another 
can  make  the  calculation  of  the  matrix  components  or  the  solving  of  the 
secular  equation  more  convenient,  but  the  final  result  is  the  same  regard¬ 
less  of  the  set  of  atomic  orbitals  used  to  start  with. 


So  far,  the  general  outline  of  the  LCAO  method  has  been  described. 
When  treated  as  an  interpolation  tool,  the  complicated  rigorous  applica¬ 
tion  of  the  method  can  be  simplified.  The  first  simplifying  step  involves 

the  atomic  orbitals.  The  Bloch  sums  of  the  atomic  orbitals  <f>  on  the 

n 

varous  atoms  of  a  unit  cell  are  not  orthogonal  to  each  other.  This  is 
because  the  atomic  orbitals  of  different  atoms  are  not  orthogonal.  This 
complicated  situation  is  removed  by  obtaining  a  new  set  of  atomic  orbitals 
which  are  linear  combinations  of  the  original  orbitals  and  are  orthogonal 
to  each  other.  The  orthogonal  atomic  orbitals  are  denoted  by  I y  .  This 
simplification  does  not  solve  the  periodic  potential  problem.  The  non¬ 
diagonal  matrix  components  of  energy  between  the  Bloch  sums  from  differ¬ 
ent  1 b  remain  to  be  found.  The  1 1>  ,  also  known  as  Lowdin  functions, 
n  n 

have  the  symmetry  properties  of  the^  from  which  they  were  derived,  which 
thus  simplifies  the  problem  of  solving  the  secular  equation  (Ref  52:1524). 

Bloch  sums  are  constructed  from  the  Lowdin  functions  in  the  form 
(Ref  52:1501): 

If1  £®i>  exp  i(k*R±)  ^n(r-R±)  (82) 

which  are  normalized  and  orthogonal.  Periodic  boundary  conditions  are 
assumed  and  N  equals  the  number  of  unit  cells  in  the  crystal.  The  sum 
is  over  the  N  unit  cells  and  the  normalization  is  over  the  repeating 
crystal.  The  Bloch  function  is  formed  from  just  one  atomic  orbital  per 
unit  cell  which  will  have  implications  in  the  diamond  structure.  The 
next  step  is  to  find  the  matrix  component  of  energy  between  two  of  these 
Bloch  sums.  The  matrix  component  is  (Ref  52:1501): 


N-1  exp  ik»(Ru-Ri) 


j 


)H  0  (r-R.)dv 
m  j 


(83) 


where  H  the  Hamiltonian  operator,  R^  ranges  over  the  positions  of  the 

atoms  on  which  the  orbitals  0^  are  located,  and  R^  ranges  over  the 

positions  of  orbitals  0  The  sum  is  over  N.  The  double  sum  can  be 

m 

written  as  a  single  summation  because  the  single  summation  is  multiplied 
by  N,  cancelling  N  If  R^  represents  the  position  of  the  atom  in  the 
central  unit  cell  having  orbital  0^,  and  the  sum  is  taken  over  R^,  the 
remaining  single  sum  over  R^.  amounts  to  summing  over  all  the  neighbors 
of  the  original  atom.  This  gives  (Ref  52:1501): 


^(Rjexp  ik*(Rj-R±)  ^0*0 r-R.)H  ^(r-R.)dv 


(84) 


where  R^-R^  represents  the  vector  displacement  from  the  atom  on  which  the 
orbital  0^  is  located,  to  a  neighboring  atom  on  which  the  orbital  0^ 
is  located.  Each  term  depends  on  the  propogation  constant  k  and  each 
term  of  the  summation  is  related  to  a  pair  of  orbitals  of  neighboring 
atoms. 

Finding  the  Lowdin  functions  from  the  atomic  orbitals  can  be  very 
difficult;  for  the  moment,  it  will  be  assumed  to  have  been  accomplished. 
Since  each  0^  is  a  combination  of  atomic  orbitals  on  many  nearby  atoms, 
the  integral  in  Equation  (84)  represents  a  linear  combination  of  many 
Integrals  of  the  form  (Ref  52:1501): 

* 


/ 


0n  (r-R1)H^(r-R1)dv. 


(85) 
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The  Hamiltonian,  H  involves  a  periodic  potential  which  can  be  written  as 
the  sum  of  a  kinetic  energy  operator  and  a  potential.  The  potential  is 
approximately  a  sum  of  spherically  potential  wells  located  at  all  the 
atoms  of  the  crystal.  The  integral  can  then  be  viewed  as  a  linear  comb¬ 
ination  of  integrals  of  a  product  of  three  functions  (Ref  52:1501): 

. _ .  ^ — 

1)  an  atomic  function  <t>  (r-R. )  located  on  the  atom  at  R . , 

n  i  1 

2)  another  atomic  function  <t>  (r-RT")  on  the  atom  at  R.,  and 

ml  J 

3)  a  spherical  potential  function  located  on  a  third  atom. 

These  three-center  integrals  are  very  difficult  calculations  as  is  the 
calculation  of  the  matrix  components  of  energy;  in  fact,  the  calculations 
are  impossible  without  the  use  of  high-speed  computers.  This  is  the 
point  at  which  the  simplifications  of  the  LCAO  are  introduced  which  con¬ 
verts  the  method  from  a  rigorous  one  to  a  simple  interpolation  method. 

The  form  of  Equation  (84)  is  used  for  the  matrix  components  of  energy 
but  the  integrals  are  replaced  by  disposable  constants  chosen  to  fit 
accurately  determined  energies  at  particular  values  of  k.  Some  care  must 
be  taken  regarding  the  precision  of  the  method  which  depends  upon  the 
finite  number  of  disposable  constants.  Only  a  finite  number  of  arbitrary 
constants  is  necessary  to  fit  a  finite  number  of  accurate  calculations  of 
energy.  In  choosing  a  finite  number  three  simplifying  observations  are 
made  (Ref  52:1501-1504)  and  are  discussed  in  the  following  paragraphs. 

The  integrals  in  Equation  (84)  will  get  smaller  numerically  as  the 
i-th  and  j-th  atoms  get  farther  apart.  If  the  functions  iff  were  atomic 
orbitals  then  the  integral  would  go  to  zero  unless  the  two  atoms  were 
close  enough  in  proximity  to  have  their  atomic  orbitals  overlap  apprec¬ 
iably.  Since  the  are  actually  orthogonalized  Lowdin  functions 
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containing  contributions  from  the  orbitals  of  nearby  atoms,  the  integrals 
will  be  appreciable  for  greater  interatomic  distances  than  if  the  inte¬ 
grals  were  formed  directly  from  the  atomic  orbitals.  All  the  same,  the 
integrals  will  decrease  fairly  rapidly  as  the  distance  between  atoms 
decreases  (Ref  52:1501-1502). 

The  symmetry  properties  of  the  crystal  suggest  that  the  integrals 
in  Equation  (84)  are  related  (Ref  52:1502).  For  example,  symmetric  ro¬ 
tation  and  inversion  operations  about  an  atom  at  R^,  where  the  orbital 

0  is  located,  will  carry  the  atom  at  R,  into  another  similar  atom  into 
«  3 

the  same  unit  cell  or  another  unit  cell.  The  relations  between  integrals 

involving  orbitals  Ip  and  p  on  various  atoms  will  be  all  the  same  at 

n  m 

equal  distances  from  the  first  atom.  Noting  this,  the  number  of  indepen¬ 
dent  integrals  decreases  and  in  some  cases  the  finite  number  of  disposable 
constants  remaining  are  nearly  the  right  amount  needed  to  fit  the  results 
of  accurate  calculations  at  the  symmetry  points  in  the  Brillouin  zone 
which  have  already  been  made.  Using  these  symmetry  properties  results 
in  simple  formulas  for  the  matrix  components  of  energy  (Ref  42) . 

A  second  simplification  is  made  by  considering  only  the  atomic  or¬ 
bitals  whose  energy  is  somewhere  near  the  energy  of  the  energy  bands  of 
interest  (Ref  52:1502),  For  example,  when  considering  the  valence  and 
conduction  bands  of  diamond,  for  an  element  in  the  3d  transition  group, 
one  can  consider  the  3d,  4s,  and  4p  atomic  orbitals  and  disregard  the 
others.  Using  only  a  few  atomic  orbitals  reduces  the  order  of  the  sec¬ 
ular  equation  to  one  that  can  be  handled  more  easily.  In  diamond,  there 
are  eight  orbitals:  the  2s  orbital  and  the  threefold  degenerate  2p  orb¬ 
itals  for  each  of  the  two  atoms  in  the  unit  cell,  which  produces  an  eighth- 
order  secular  equation. 
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The  third  simplification  is  called  the  two-center  approximation 
(Ref  52:1502-1504).  Three-center  integrals  are  smaller  than  tws-center 
integrals  and  in  some  cases  the  former  can  be  disregarded  thus  simpli¬ 
fying  the  equation  even  further.  When  the  potential  energy  in  the  Ham¬ 
iltonian  is  written  as  the  sum  of  spherical  potentials  located  on  the 
various  atoms,  the  three-center  integrals  can  be  neglected.  Part  of  the 
potential  energy  remains  which  represents  the  sum  of  spherical  potentials 
located  on  the  two  atoms  on  which  the  atomic  orbitals  are  located,  anal¬ 
ogous  to  the  case  of  a  diatomic  molecule.  Let  R_.-R^  represent  the  vector 
from  one  atom  to  another,  i.e.  the  axis.  The  functions  Ip  can  be  expanded 
as  linear  combinations  of  the  components  of  angular  momentum  around  the 
axis:  Of  7T,  and  5.  Although  the  (p's  are  Lowdin  orthogonalized  functions 
and  not  atomic  orbitals,  they  have  the  same  symmetry  properties  with  res¬ 
pect  to  the  crystal  so  the  expansion  is  still  a  good  approximation. 

Each  integral  can  be  symbolized  by  ,m,n,) ,  for  example,  where  l,m 

and  n  denote  the  direction  cosines  of  the  R^-R^  vector  (Ref  52:1503). 

The  first  subscript  x,y,  or  z  refers  to  the  p^,  p^,  or  p ^  functions, 
respectively.  The  second  subscript  refers  to  the  various  d  functions; 

these  subscripts  stand  for  the  functions  whose  dependence  on  angle  is  of 

2  2  2  2 

the  form  xy,  yz,  xz,  x  -y  ,  or  3z  -r  multiplied  by  the  appropriate 
function  of  r.  The  energy  integral  represents  an  integral  in  which  the 
function  (p  is  a  p  -like  function,  and  lp^  is  a  d  function  with  symmetry 
properties  like  xy.  Ex»  (1 ,m,n, )  is  a  function  which  can  be  written  in 
terms  of  two  integrals:  the  first  is  between  a  p <3  orbital  on  the  first 
atom  and  a  d<7  orbital  on  the  second,  symbolized  by  (pd<7);  and  the  second 
is  between  a  p7T  orbital  on  the  first  atom  and  a  d 7T  orbital  on  the 
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TABLE  II 


Energy  integrals  for  a  crystal  in  terms  of  two  center  integrals 


E,., 
E,., 
£*. 
Et.it 
Et.  i’-i* 
£n.’-r* 

Eg.  n 

Et., i 

E,.  tl 


fi.ii'-f1 

£o.*-r« 

A«». » 

E,,iK,> 
Eu.iK,* 
E,,.  liLr* 


(ssa) 

l(spo) 

P(ppa)+(1—P)(ppx) 

Imippa)  —Im(ppx) 
ln(ppa)-ln(ppx) 

\3  Im(sda) 

\dl(T—nP)  (Jiiff) 

+  —  2P)(pdr) 

\llmn!  pda)  —  Hmn(pdx) 

\3Pn(pda)+n(\-2P)(pdx) 

\\V(P- m*)  ( pda)  +  l(  1  -P+m*)(  />rfr) 

J  v3  m  (P — m* )  ( pda }  —  m  ( 1+  P  —  m5)  ( pdr ) 
i\3n(P-m,)( pda)  -„tp -„'■)( pdT) 

—  1  (f*+  «p)  —  v3/nJ  f  pdr) 

(n*  —  1  (^4-  w1)  3  f  pda)  —  v3rrr.J ;  pdx) 

«C«5  —  }  fP+  op)  2(pda) 4- \3»  (/*+  ms)  ( /></*) 
iPjrd  (dda)  -f-  (p+m-  —  iPmi)(ddr)  4-  (n-+PnP)(ddP) 

Hnpn[dda)  -M»(l  —  4w2)  (ddx)Jrln(np  —  1  ''(dd&) 

3Pmn(dda)-rmn(l  —4P)(ddx)P-mn(F—\)(ddS) 

>,lm(P- m'-)  (dda) +2  lm(m'—P^  (ddx) +\lm(P-  m!)  (ddt) 
jmn(P-m')(dd*)-mr.[l  +  2(F-m>)2(ddr)+mnlU-\(P-m')2(ddS) 
lnl(P-m:)(dda)-i-r.lll-2(P-m=)2(<!dr)-nltl-}(P-m')yddi) 

'JSlm[n*  —  J  (P  -f  m!)  ]  (&/«•)  —  2  v3  ImtP  (ddw) 4-  J'3/m  ( 1 4-  n*)  (ddt) 

vSmfifn1  —  J(P4-  m')’}(dda)+\'3mn(F+m*  —  n:)(ddr)  —  ]\‘$mn(P+m*)(ddi) 

\Vntn’-i(P+m’)2(dda)+\yn(P+m'-n')(ddT)-Wln(P+m>)(ddi) 

l  (P  -  m’)’  (dda) + [P+ m*  -  (P  -  m2)*]  (ddx) + [n=+  }  (P  -  m2)«]  (ddi) 

\'4(P-m’)[n'-l(P+mi)-](dda)+\l3,P(m’~F)(ddx)+{\'3(l+n’-)(P-m*)(ddi) 

[»5  -}(/»+  m’)j>(dda)  4-  3«s  (P+ m»)  (ddx)  +l(P+mrlt(ddi) 


(Ref  52:1503). 

second  atom,  symbolized  by  (pd7T).  The  first  index  refers  to  the  first 
orbital  and  similarly  for  the  second.  The  energy  integral  was  determined 
by  the  product  of  the  three  integrals  outlined  earlier.  Formulas  for  the 
energy  integrals  for  all  the  combinations  of  the  s,  p,  and  d  functions 
are  given  in  Table  II  (Ref  52:1503).  By  using  Table  II,  the  various 
energy  integrals  E  can  be  approximated  in  terms  of  a  smaller  number  of 
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integrals  which,  is  a  good  approximation  in  some  cases  (Ref  52:1504).  The 
table  is  also  very  useful  in  making  comparisons  with  other  existing  calc¬ 
ulation  which  automatically  assume  two-center  integrals. 

Applying  the  LCAO  to  the  Diamond  Structure  (Ref  52 : 1516-1522)  .  The 
diamond  structure  can  be  described  as  a  face-centered  cubic  structure 
having  two  atoms  in  the  unit  cell  rather  than  one.  One  atom  is  located 
at  pai  +  qaj  +  rak,  site  1,  and  the  other  atoms  are  displaced  by 
h  a(i+j+k),  called  site  2,  where  p,  q,  and  r  are  integers,  p+q+r  is  an 
even  integer,  and  a  is  the  lattice  constant.  Eight  Bloch  sums  are  ob¬ 
tained  by  setting  up  Bloch  functions  of  s,  p  ,  p  ,  orbitals  on  each 
of  the  two  types  of  atomic  sites.  The  matrix  components  of  energy  obtained 
by  means  of  the  Bloch  sums  are  denoted  by  symbols  such  as  (s/x)^*  This 
represents  a  matrix  component  between  an  s-orbital  Bloch  function  on  site 
1  and  a  p^  Bloch  function  on  site  2.  There  are  four  nearest  neighbors 
of  an  atom  in  the  diamond  structure  as  shown  in  "igure  21  (Ref  35:9). 

If  the  atom  under  consideration  is  at  the  origin,  its  neighbors  are  lo¬ 
cated  at  the  positions  a(%  h  h)  ,  a(-b  ~h  ~h)  >  a(~h  %  ~h) ,  and  a(-%  -h  h) . 
The  second  nearest  neighbors  are  the  type  denoted  by  a(110).  The  resul¬ 
ting  matrix  components  are  given  in  Table  III  (Ref  52:1517)  where  £,  T\ , 
and  £  denote  the  symmetry  sumbols  for  the  diamond  structure  and  the 
energy  Integrals  (Ref  52:1517)  are  given  by 

£*.(iiJ)=  («»)». 

£*,(ili)  =  J  (PP<r)  I- \  (PP*)  1. 

£«..(*}!)= 3"1  (*^)i, 

£.,«(110)  =  (sse)t,  (86) 

£*.(011)=  (pfirh 
E,.,(m  =  l(PP*h-l(PP'h 

Et,(110)  =  }(pp<T)t+l(PPr)*> 

£*.(uo)=r»(^)». 

£»  .(01 1)  =  £«, ,  (01 1)  =  0. 
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Fig  21.  The  diamond  lattice  (Ref  29:28;  35:9). 

The  resulting  Hamiltonian  matrix  is  real  and  symmetric,  and  has  the 
following  components  (Ref  52:1517): 


(j/j)u=  (*A)»» 

(z/x)n=  (*/*)« 

(j/i)u=  (»A)»i* 

(j/*)i»“  -  (*/*)ii* 

(s/x)it-  —  (*A)«* 

(.r/.r)u=  (*/*)„• 

(x/y)n=  (x/y)it*=  (yA).j 

f.t/y)n=  (*/)•)«* 


,  (000) -f  .(110)  (cosf  cosi+cosij  cosf +cosf  cosf) 

. (000)  +  4£., . (1 1 0) (cosf  cost;!- cosf  cosf)+4£.,«(011)  cost;  cosf 
4-En(jjj)(cosjf  cosJt;  cosjf — »  sinjf  sinh  sinjf) 
iE,;,(j  j  j)(t  sinjf  cosh  cosjf— cosjf  sinh  sinjf) 

—  4£,.,(011)  sinf  sini;+4:£,,I(110)(sinf  cosn+sinf  cosf)  (87) 

4.E.,  •  (Jjj)  (cosjf  cos  jv  cosjf sinjf  sinh  sinjf) 

4E».n(j  j  j)(«  cosjf  cosjTj  sin  jf— sin  jf  sinji;  cosjf) 

—  4£.,,(110)  sinf  sini;  — 4i£.,,(01I)(sinf  cosf — sim;  cosf) 


This  can  be  reduced  to  a  2x2  matrix  whose  determinant 


Bll+H13-E 


H12+H14 


H12+H14 


H22+H24-E 


(88) 


assuming  A  symmetry  (Ref  52:1517). 
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For  the  (100)  and  (111)  directions,  the  equation  predicts  a  number  of 
levels  and  their  degeneracy,  and  the  symmetry  at  the  edges  of  the 
Brillouin  zone.  Some  values  of  the  energy  integrals  vanish  as  a  result 
of  the  two-center  approximation  but  the  same  type  of  integrals  can  be 
chosen  to  approximate  known  results.  Much  use  has  been  made  of  the 
diamond  structure,  partly  due  to  Herman's  results  using  the  OPW  method 
(Ref  24;  25)  and  the  comparison  with  the  results  of  the  LCAO  method. 

The  general  agreement  comes  from  the  minimal  contributions  of  the  second- 
nearest  neighbors  compared  to  those  of  the  nearest  neighbors.  The  OPW 
method  was  applied  to  diamond  by  Herman  (Re  24),  to  Ge  by  Herman  and 
Callaway  (Ref  8:74-78;  23),  and  to  Si  by  Woodruf  and  Bassani  (Ref  8:70-74). 
The  results  are  shown  in  Figures  22  and  23. 

Suggestions  for  Applying  the  LCAO  Method  to  the  Ge-Si  Crystal . 

There  are  many  things  to  consider  in  applying  the  LCAO  method  to  the 
Ge-Si  alloy  of  interest.  The  alloy  is  not  a  pure  crystal  by  definition. 

In  addition  to  the  binary  property  of  the  alloy,  the  microscopically 
nonuniform  distribution  of  each  element  does  not  facilitate  easy 
application  of  the  LCAO  method.  Assuming  the  lattice  retains  its 
crystal  proerties  at  10%  Si  and  90%  Ge,  the  uneven  nature  of  the  comp¬ 
osition  suggests  very  complicated  Bloch  sums  and  thus  very  complicated 
wavefrnctions  The  binary  alloy  keeps  its  diamond  crystal  structure 

but  a  difficulty  arrises  in  representing  a  10%  Si  and  90%  Ge  alloy  in  a 
tetrahedral  unit  cell. 

A  recommended  approach  is  to  take  the  statistical  average  of  the 
results  of  several  applications  of  the  LCAO  method  to  a  crystal  which  may 
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Fig  22.  a)  The  energy  bands  of  a  diamond  crystal  in  the  100  and  111 
directions  (Ref  24) , 

b)  the  energy  band  structure  of  a  Ge  crystal  (Ref  25), 

c)  the  energy  bands  of  Ge  along  the  100  and  111  axes  (Ref  8:76) 
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Fig  23.  The  energy  bands  of  Si  along  the  100  and  111  directions  (Ref  8:72). 

.  be  10%  Si  over  its  total  volume  but  not  over  the  volume  of  the  unit  cell. 

The  LCAO  method  involves  linear  combinations  of  Bloch  sums  between 
atomic  orbitals  of  atoms  even  though  the  atoms  under  consideration  may 
not  be  the  same.  Assume  also  that  the  periodic  potential  remains 
periodic  after  the  substitutional  addition  of  Si  in  the  Ge  crystal  is 
10%  Si  and  90%  Ge,  the  distribution  is  a  uniform  as  possible,  and  the 
location  of  the  Si  atoms  are  periodic  when  moving  from  one  unit  cell  to 
another  (  a  reasonable  assumption  for  crystals  grown  by  epitaxy) ,  then 
the  LCAO  method  can  be  applied  to  each  unit  cell,  producing  an  energy 
band  structure  for  each  cell  in  the  volume. 

It  is  obvious  that  a  single  unit  cell  cannot  have  a  10%  composition 
of  Si  because  whole  atoms  must  be  preserved.  The  binary  alloy,  then,  is 

( 
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really  made  up  of  unit  cells  which  may  have  zero,  one,  two,  or  more  Si 
atoms  mixed  with  the  Ge  atoms  of  the  unit  cell  so  that  the  concentration 
over  the  volume  of  the  crystal  will  have  10%  Si  -  90%  Ge.  When  the 
LCAO  method  is  applied  to  a  unit  cell  having  no  Si  atoms,  the  resulting 
energy  hand  structure,  which  has  already  been  calculated  by  Herman  and 
Slater  and  Koster  (Ref  25;  52),  will  be  that  of  Ge  as  shown  in  Figure  22. 
When  the  LCAO  method  is  applied  to  unit  cells  having  one  Si  atom,  the 
resulting  energy  band  structure  will  closely  resemble  that  of  Ge  with 
slight  modifications  of  the  Si  (100)  and  Ge  (100)  conduction  band  curves 
(Ref  5;  8:76;  27).  Similar  applications  of  the  LCAO  mathod  would  have 
to  be  done  for  all  the  unit  cells  of  the  crystal  or  a  representative 
subgroup.  A  weighted  average  of  the  resulting  energy  band  structures 
should  produce  an  energy  band  structure  which  describes  the  macroscopic 
behavior  of  holes  and  electrons  in  the  Ge-Si  alloy. 

It  would  be  helpful  in  evaluating  the  constants  arising  in  the  LCAO 
results  to  have  some  idea  of  the  energy  band  structure  along  the  symmetry 
lines  of  the  crystal.  This  information  can  be  obtained  from  the  results 
of  optical  and  cyclotron  resonance  experiments  (Ref  35:290-307;  29:318- 
322;  59)  as  in  the  results  of  Johnson  and  Christian  for  a  0-199%  range 
of  Si  in  Ge  (Ref  27).  The  absorption  properties  of  the  alloy  have  been 
studied  by  Braunstein,  Moore,  and  Herman  (Ref  5)  for  the  same  compositions 
of  Si  and  Ge.  At  about  15%  Si  the  energy  band  structure  changes  from  a 
Si-like  (100)  to  a  Ge-like  (100)  conduction  band  structure  which  is  the 
crucial  point  of  interest  for  the  subject  of  this  thesis. 

( 
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While  work  progresses  in  the  area  of  the  crystal  growth  of  the 
Ge~Si  alloy,  the  predicted  results  of  the  quantum  mechanical  analysis 
of  the  binary  alloy  forms  the  basis  for  the  decision  to  treat  the 
photodetector  alloy  as  pure  Ge  with  regard  to  its  chemical  and  material 
properties  which  come  into  play  during  the  fabrication  process.  The 
following  chapter  describes  the  fabrication  procedure  for  the  five  Ge 
photodiodes  discussed  in  the  thesis. 
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IV  Device  Fabrication  and  Testing 

The  five  Ge  photodiode  samples  which  were  analyzed  in  this  thesis 
research  are  secribed  in  this  section  in  terms  of  device  fabrication. 

The  process  for  each  sample  can  be  found  in  Table  III.  The  three  samples 
that  were  implanted  are  G-7-1,  G-7-4,  and  G-8-1 .  The  diffused  samples, 
D-6-1  and  D-6-3,  were  separate  diodes  from  the  same  wafer.  The  starting 
material  for  all  samples  was  3  f2-cm,  p-type  Ge  wafers  three  inches  in 
diameter. 

The  implanted  diodes  were  processed  as  130  x  135  mil/square  chips 
which  underwent  processing  mounted  in  an  implant  holder  shown  in  Figure 
24.  The  holder  was  capable  of  processing  four  chips  each  containing 
four  diodes  in  a  square  matrix.  The  30-mil  implant  mask  is  also  shown 
in  Figure  24  next  to  the  3-mil  contact  mask  fitted  inside  the  implant 
holder. 

In  the  case  of  implanted  samples,  the  individual  chips  were  sputtered 

o 

with  Si^N^  for  five  minutes.  This  1000  A,  very  thin  film  provided  a  pro¬ 
tective  film  on  the  surface  of  the  chip  during  the  following  implant  step. 
Residual  hydrocarbons  from  the  vacuum  system  of  the  accelerator  can  be 
polymerized  on  the  chip  surface  by  ion  implantation;  this  is  analogous 
to  hitting  photoresist.  The  polymerized  hydrocarbons  can  interfere  with 
the  implantation  and  are  also  very  hard  to  remove.  The  As  implants  were 
performed  as  listed  in  the  process  schedule  of  Table  III,  followed  by 
thermal  annealing. 

Front  contacts  of  Ni/Au/Ge  were  evaporated  on  the  chips  using  a 
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3-mil  lift-off  mask  to  define  the  contact  pattern.  The  procedure  was  as 
follows.  A  mask  containing  an  array  of  3-mil  holes,  as  shown  in  the  im¬ 
plant  holder  of  Figure  24,  was  used  to  cut  holes  in  the  photoresist  as 
shown  in  Figure  25(a).  A  commercial  buffered  HF  etch  was  used  to  remove 

the  Si^N^  inside  the  3  mil  holes  in  Figure  25(b).  After  the  Si^N^  was 

o 

lifted  off  the  wafer  inside  the  3  mil  holes,  1000-1200  A  of  Ni/Au/Ge 
metallization  was  evaporated  on  the  wafer  as  shown  in  Figure  25(c).  The 
thick  photoresist  (AZ)  layer  enabled  a  discontinuity  to  form  between  the 
metal  that  forms  the  contact  and  the  metal  deposited  on  top  of  the  photo¬ 
resist.  The  discontinuity  allowed  the  latter  to  be  lifted  off  at  the 
same  time  the  photoresist  was  removed  with  acetone.  The  remaining  struc¬ 
ture  is  shown  in  Figure  25(d)  with  a  final  metallization  thickness  of 
o 

1000-1200  A. 

After  back  contacts  were  evaporated  on  the  chips,  two  samples,  G-7-1 
and  G-8-1 ,  were  left  planar  and  G-7-4  was  mesa  etched.  Dots  of  black  wax 
were  placed  by  hand  over  the  contacts  of  a  sample  and  part  of  implanted 
region  as  shown  in  Figure  26(a) .  The  sample  was  then  mounted  on  a  glass 
slide  and  placed  in  a  peroxide  etch  (1  :  5  DI-H^O)  for  43  minutes. 

The  peroxide  etches  at  a  rate  of  0.0275  }lm  per  minute  (Ref  46:88-91).  A 
rinse  in  acetone  following  the  etch  removed  the  wax  leaving  a  mesa  struc¬ 
ture  as  shown  in  Figure  26(b).  Obviously,  diode  area  by  this  means  was 
not  well  controlled.  The  result  of  the  imprecision  will  be  discussed  in 
Chapter  VII. 

A  number  of  diodes  were  probed  using  the  set-up  in  Figure  27  where 
the  probes  were  connected  to  a  conventional  curve-tracer.  Chips  containing 
working  diodes  were  mounted  on  a  28-pin  DIP  package.  Ultrasonic  bonding 
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Fig  26.  Mesa  etching  using  dots  of  black  wax  over  the  diode  area. 

at  room  temperature  damaged  some  samples  which  explains  the  difference 
in  the  bonding  of  sample  G-8-1.  One  problem  with  the  shallow  junction 
is  that  the  "shallowness"  is  very  sensitive  to  the  pressure  placed  on  it 
by  bonding  wires  and  bonding  procedures  themselves.  Some  samples  that 
displayed  excellent  I-V  curves  prior  to  bonding  were  unfortunately  lost 
as  a  result  of  being  shorted  by  ultrasonic  bonding.  By  heating  the  pack¬ 
age  to  280°C,  it  was  possible  to  lower  the  power  and  time  settings  for 
ultrasonic  bonding.  For  the  four-diode  chips,  bonds  were  made  from  the 
contacts  to  pins  #6,  9,  20,  and  23.  A  ground  wire  was  also  bonded  from 
the  substrate  to  pin  #15.  See  Figure  28. 

The  diffused  samples  D-6-1  and  D-6-3  originated  from  the  same  3  inch 
wafer.  The  3  ft-cm  p-type  Ge  wafer  experienced  diffusion  of  phosphorus 
from  a  phosphene  ambient  for  60  minutes  at  775°C.  Front  contacts  were 
evaporated  and  patterned  using  a  conventional  mask  consisting  of  an  array 
of  3-mil  holes.  The  lift-off  process  described  in  Figure  25  was  used 
here  also.  Back  contacts  were  omitted  in  favor  of  mounting  chips  on  the 
DIP  using  silver  paint  because  of  difficulties  with  the  vacuum  system  in 


15  ground. 


the  evaporator.  A  mask  containing  an  array  of  20-mil  holes  was  aligned 
with  the  3-mil  contacts  to  pattern  the  phtoresist  for  mesa  etching.  The 
remaining  structure  was  similar  to  that  shown  in  Figure  26.  After  scribing, 
only  one  diode  has  contained  in  each  chip.  Bonding  wires  were  connected 
to  pins  6  and  15  for  the  contact  and  ground  respectively.  A  close-up  view 
of  a  single  diode  is  shown  in  Figure  29.  It  shows  the  center  contact, 
the  bonding  wire,  and  the  perimeter  of  the  diode  defined  by  the  mesa  etch. 

The  diffusion  process  followed  the  same  sequence  as  the  implantation 
schedule  except  that  the  diffused  diodes  were  processed  as  complete  wafers 
and  the  P  impurity  atoms  for  the  n+  region  were  supplied  by  the  phosphene 
ambient  in  the  diffusion  furnace.  In  contrast  to  the  implanted  samples, 
mesa  etching  was  required  in  the  diffused  samples  to  define  the  perimeter 
of  the  area  of  the  diode.  The  mesa  etch  was  dome  using  photoresist  and 
masking  steps  instead  of  the  dots  of  black  wax  described  earlier.  Pack¬ 
aging  and  testing  of  the  diffused  samples  was  the  same  as  for  the  implant¬ 
ed  samples. 

Table  III  summarizes  the  device  fabrication  steps  of  the  implanted 
and  diffused  Ge  photodiodes.  The  finished  samples  were  probed  and  pack¬ 
aged,  then  tested  again  for  reassurance  that  workable  diodes  still  exist¬ 
ed  after  bonding.  It  should  be  noted  that  the  performance  of  some  of  the 
Ge  photodetectors  as  diodes  degraded  somewhat  after  bonding.  The  next 
step  was  to  test  the  performance  of  the  Ge  photodetectors  for  their  opti¬ 
cal  characteristics,  namely  photoresponse.  This  requires  a  discussion  of 
the  photometric  system  used  to  test  the  detectors  which  follows  in  the 
next  chapter. 


81 


V  Experimental  Procedure  and  Equipment 


Photometric  System 

The  photometric  system  used  to  test  the  photoresponse  of  the  samples 
is  shown  in  block  diagram  form  in  Figure  30  (Ref  6:3).  It  requires  bas¬ 
ically  four  components:  source,  detector,  amplifier,  and  recorder.  The 
source  emits  optical  radiation  over  a  wide  range  of  wavelengths,  which 
is  incident  on  the  photodetector.  The  photodetector  does  not  respond 
over  the  same  wide  spectrum  of  wavelengths  emitted  by  the  light  source. 

In  the  smaller  range  of  wavelength  over  which  it  does  respond,  incident 
light  generates  enough  electron-hole  pairs  to  produce  a  very  weak  signal. 
The  weak  output  signal  is  amplified  to  appropriate  levels  for  signal  proc¬ 
essing  and  is  recorded. 

Weak  Signals  Buried  in  Noise.  There  are  various  noise  limitations 
associated  with  the  detection  of  optical  radiation.  The  main  problem  is 
that  the  eceltrical  signal  is  so  weak  it  is  usually  buried  in  noise.  One 
requirement  to  obtain  meaningful  measurements  is  to  reduce  the  system 
noise  to  that  of  the  detector  or  the  source.  The  detector  must  have  at 
least  three  times  the  noise  resistance  of  the  amplifier  (Ref  6:10).  A 
comparison  of  three  systems  is  made  in  Reference  6  and  the  synchronous 
detection  scheme  described  in  Appendix  A  (Ref  6;  26;  45;  50)  is  preferred 
over  the  DC  and  conventional  AC  systems.  The  preference  can  be  attributed 
to  minimized  noise  problems. 

Equipment .  In  terms  of  the  block  diagram  of  the  photometric  system 
shown  in  Figure  30,  the  radiant  source  takes  the  form  of  a  globar.  The 
globar  is  a  water-cooled,  silicon  carbide  light  source  whose  spectral 
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Fig  30.  The  block  diagram  of  a  photometric  system  (Ref  6:3). 


characteristics  have  been  calibrated  and  measured  for  operation  at  230 

watts  of  globar  power.  The  characteristics  are  listed  in  Table  IV  by 

drum  number,  the  corresponding  wavelength  A(//m),  energy  E(eV) ,  and 

2 

globar  power  per  unit  area  H(^W/cm  ).  The  drum  number  comes  from  the 
monochromator  used  to  manually  tune  the  desired  wavelength  or  scan  the 
spectrum  of  the  globar  light  source. 

Erratic  behavior  on  the  part  of  the  globar  light  source  prompted  the 
run  of  a  second  set  of  calibration  data.  The  first  set  of  data  was  taken 
on  May  30,  1979.  During  the  course  of  the  thesis  work,  the  globar  was 
stable.  Instability  in  the  operation  of  the  globar  near  the  end  of  the 
thesis  quarter  required  the  recalibration  on  September  28,  1979;  the  un¬ 
stable  behavior  was  due  to  contact  problems.  The  second  set  of  calibra¬ 
tion  data  tended  to  flatten  the  output  power  spectrum  of  the  globar  in  a 
nonlinear  fashion  as  can  be  seen  by  the  column  of  Table  IV. 


TABLE  IV 


Calibration  Data 


Drum 

Number 

E 

(eV) 

A 

(fim) 

H* 

(^W/cm^) 

H** 

(JM/cm  ) 

AH  2 
( flV/cm  ) 

1580 

1.19 

1.04 

28 

22 

6 

1570 

1.11 

1.11 

140 

142 

2 

1560 

1.04 

1.19 

307 

306 

1 

1550 

0.973 

1.28 

432 

437 

5 

1540 

0.908 

1.37 

523 

516 

7 

1530 

0.852 

1.46 

643 

623 

20 

1520 

0.800 

1.55 

721 

695 

26 

1510 

0.752 

1.65 

766 

730 

36 

1500 

0.710 

1.75 

777 

742 

35 

1490 

0.672 

1.85 

750 

713 

43 

1470 

0.611 

2.03 

768 

728 

50 

*  Calibration 

data  taken 

30  May 

1979. 

**  Calibration  data  taken  28  September  1979. 

The  monochromator  is  a  Model  098  (Serial  No.  297)  from  the  Perkin- 
Elmer  Corporation  in  Norwalk,  Ct.  It  is  a  single  pass,  basic  prism  mono¬ 
chromator  designed  for  use  in  the  infrared  region  but  can  be  used  with 
other  prisms  to  change  the  wavelength  region  or  vary  the  amount  of  dis¬ 
persion.  Wavelength  is  scanned  manually  by  rotating  a  drum  which  adjusts 
the  internal  optics  of  the  monochromator  to  vary  the  wavelength.  A  sche¬ 
matic  diagram  of  a  monochromator  is  given  in  Figure  31. 

Figure  32  shows  the  experimental  set-up  for  the  photometric  system. 
The  tall  cylindrical  shaft  is  the  vacuum  chamber,  the  top  of  which  holds 
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WINDOW 


monochromator  (Ref  39) . 


Fig  32.  The  experimental  set-up  for  the  photometric  system. 

the  DIP  package  on  which  the  photodiode  chip  is  mounted.  Refer  to  the 
diagram  in  Figure  33.  The  diode  is  clamped  to  the  header  and  the  cylinder 
slides  over  the  header  making  a  seal  for  pumping  down  the  chamber.  The 
vacuum  gauge  at  the  top  center  of  the  photo  in  Figure  32  monitors  the 
vacuum  pressure  inside  the  chamber.  Coaxial  cables  from  the  photodetector 
lead  to  a  preamplifier  and  to  the  rack  of  instruments  shown  in  Figure  34. 
The  instruments  are  listed  in  Table  V  and  their  uses  will  be  demonstrated 
in  the  following  paragraph. 

System  operation  can  be  summarized  as  follows.  The  monochromator 
selects  a  wavelength  provided  by  the  globar  and  emits  the  infrared 
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Fig  33.  Laboratory  set-up  for  mounting  a  sample  for  testing. 


Fig  34.  Panel  of  instruments  used  in  processing  and  recording  the  output 
of  the  photodetector. 

radiation  through  a  slit.  The  sample  mounted  on  the  vacuum  header  receives 
the  incident  light  through  a  silicon  glass  window  which  is  aligned  with 
the  output  window  of  the  monochromator.  The  x-y  alignment  of  the  diode 
is  adjustable  at  the  base  of  the  vacuum  chamber.  The  bias  voltage  is 
set  using  the  KEPCO  programmer  and  DC  regulated  supply.  Coaxial  cables 
from  the  diode  under  test  lead  to  a  preamplifier  which  feeds  the  lock-in 
voltmeter  for  signal  processing  as  outlined  earlier.  The  5055A  digital 
recorder  prints  the  output.  A  100K  resistor  is  connected  in  series  with 
the  diode  so  that  what  is  actually  being  measured  is  the  voltage  across 
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the  resistor  due  to  the  photocurrent  generated  at  the  wavelength  selected 
by  the  monochromator.  Therefore,  the  3490A  Multimeter  is  actually  meas¬ 
uring  output  voltage.  It  will  be  converted  to  photocurrent  later  for 
determination  of  the  photoresponse  characteristics. 

Full  details  of  the  experimental  procedure  are  described  in  Appendix 
B.  The  test  run  consisted  mainly  of  placing  the  diode  under  0,  5,  and 
10V  bias  and  obtaining  the  output  voltage  across  a  100K  resistance  for  the 
spectrum  of  wavelengths  listed  in  Table  IV.  The  raw  data  are  presented 
in  the  following  results  chapter. 


TABLE  V 


List  of 

Equipment  used 

in  the  Photometric 

System. 

Instrument 

Model 

Serial  No. 

Manufacturer 

Monochromator 

098 

297 

The  Perkin-Elmer 
Corporation 

Norwalk,  Ct. 

Oscilloscope 

RM503 

002185 

Tektronix,  Inc. 
Portland,  Oregon 

Lock-In 

Voltmeter 

131 

Brower  Laboratories 

Inc . 

Westboro,  Mass. 

Programmer 

500 

Brower  Laboratories 

Inc . 

Westboro,  Mass. 

Multimeter 

3490A 

Hewlett  Packard 

Digital 

Recorder 

5055A 

Hewlett  Packard 

ABC  Regulated 

DC  Supply 

KEPCO 

Digital 

Programmer 

KP-1 

C-64409 

KEPCO 

Digital 

Voltmeter 

3440A 

Hewlett  Packard 

Automatic  Range 
Selector 

3442A 

Curve  Tracer 

577 

B071481 

Tektronix,  Inc. 

Standard  Test 
Fixture 

177 

Tektronix,  Inc. 
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VI  Results 


This  chapter  presents  the  experimental  results  obtained  with  five 

photodetector  diodes  of  varying  structures.  The  results  include  the 

forward  and  reverse  bias  diode  characteristics,  the  diode  and  contact 

diameters,  and  the  output  voltage  as  a  function  of  applied  bias  for  each 

sample.  The  I-V  data  points  are  listed  in  Tables  C-l  through  C-5  in 

Appendix  C  by  sample  number.  The  I-V  characteristics  will  be  given  in  a 

plot  later  in  the  next  section.  A  list  of  diameters  an<.  resulting  total 

areas  is  given  in  Table  VI.  The  diameters  for  each  diode  were  obtained 

by  photographing  the  diode  at  high  magnification  (80x)  and  measured  by 

using  a  1mm  scale  photographed  at  the  same  magnification  as  shown  in 

Figure  29.  The  area  was  calculated  in  square  centimeters  for  convenience. 

For  example,  the  total  active  area  of  photodiode  G-7-1  #9  was  found  by 

subtracting  the  contact  area  from  the  total  area  of  the  diode: 

2 

total  area:  7T(0.062  cm)  /4 

2 

contact  area:  -  7T(0.Q34  cm)  /4 

-3  2 

active  area:  2.11  xlO  cm 

which  is  the  first  entry  in  Table  VI.  It  was  assumed  that  no  radiation 
penetrated  the  metallization  of  the  contact. 

As  explained  in  Chapter  V,  the  output  voltage  recorded  is  actually 
the  voltage  across  a  100K  resistor  in  series  with  the  photodetector  which 
is  due  to  the  photocurrent  generated  therin.  These  photovoltages  are 
listed  in  Tables  D-l  through  D-10  in  Appendix  D  for  each  sample  for  0,  5, 
and  10  volts  bias  for  each  set  of  calibration  data.  These  voltages  are 
used  to  determine  the  photocurrent  values  simply  by  dividing  by  the  100K 
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TABLE  VI 


Diameters  of  Diodes  and  Contacts  Used  to  Determine  Active  Area 


Total 

Diameter 

Contact 

Diameter 

Active 

Area 

G-7-1 

it  9 

0. 62mm 

0.34mm 

2.11 

x  10-3 

2 

cm 

G-7-4 

#23 

0.52ram 

0 . 33mm 

1.27 

x  10-3 

2 

cm 

G-8-1 

#23 

0 . 97mm 

0.34mm 

6.48 

x  10“3 

2 

cm 

D-6-1 

#  6 

0. 49mra 

0 . 08mm 

1.84 

x  10“3 

2 

cm 

D-6-3 

it  6 

0.49mm 

0 . 08mm 

1.84 

x  10“3 

2 

cm 

resistance.  At  5  volts,  the  output  voltage  of  G-7-1  it 9  at  a  wavelength 
of  1.46/im,  is  122225  /JV.  To  obtain  the  photocurrent,  Ohm's  law  is  applied 

1  =  V  =  1222225  UV=  1.22225  //A. 

R  100  K12 

The  values  for  photocurrent,  incident  optical  power,  responsivity  and 

quantum  efficiency  are  tabulated  in  Tables  E-l  through  E-30  in  Appendix 

E.  Since  the  active  area  varied  with  each  device,  it  was  important  to 

know  the  actual  optical  power  falling  on  a  particular  diode.  Given  this 

2 

quantity  in  fjU  per  cm  for  each  wavelength,  because  the  output  of  the 
globar  was  not  uniform  over  its  spectral  range,  it  was  a  simple  matter 
to  convert  the  optical  power  to  fi W  in  direct  relation  to  the  area  of 
a  specific  diode.  The  power  is  listed  next  to  the  photocurrent  for  ease 
in  calculating  the  photoresponse. 

As  a  sample  calculation,  consider  the  photoresponse  of  G-7-1  #9  at 

a  wavelength  of  1.46  fim  and  5  volts  bias.  It  is  obtained  as  defined  in 

2 

Equation  (7)  of  Chapter  II.  The  output  of  the  globar  is  623  H W/cm  ,  or 
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(623 MW/cm^)x(2. 111x10  ^cm^)=  1.3152 /iW  in  this  case.  Then 


R  -  I  -  1.22225  MA  =  0.929  /ttA/MW 
H  1.3152  *tW 


Quantum  efficiency  is  obtained  by  Equation  (15)  in  Chapter  II: 


7)  =  pi  =  (q  9291  — !-■  ftfrA  =  7Q  7 y 

RXCATm)  (1.46)  79,7/s 


It  is  defined  as  the  number  of  electrons  generated  per  photon. 

A  better  grasp  of  the  significance  of  the  results  of  the  experiment 
can  be  obtained  from  plots  of  photoresponse  and  quantum  efficiency  as  a 
function  of  wavelength.  These  are  shown  and  discussed  in  the  following 
section.  Other  results  which  were  previously  given  in  Table  III  but  not 
mentioned  in  context  are  discussed  in  the  following  paragraph. 

The  availability  of  spare  samples  made  it  possible  to  bevel  and  stain 
a  diffused  diode  to  determine  its  junction  depth,  x^~1.2j/m.  Although 
the  luxury  of  spare  samples  did  not  extend  to  the  implanted  samples,  an 
estimate  of  the  junction  depth  can  be  made  from  the  predictions  of 
Linhart,  Scharff  and  Schi^tt  (LSS  theory).  At  200  KeV,  an  As  implant  can 
be  expected  to  reach  a  junction  depth  of  ~  1.0  (Ref  17).  It  is 
known  from  other  work  that  for  As  implants  in  Si  and  Ge  under  conditions 
like  those  in  this  experiment  the  measured  impurity  profiles  are  reason¬ 
ably  close  to  those  of  LSS  theory.  A  measurement  of  carrier  lifetime  in 
the  diffused  diodes  came  to  T  =  2.25  p sec  which  will  also  be  discussed 
further  in  the  next  section. 

In  summary,  detailed  data  was  obtained  on  a  fairly  restricted  set  of 
samples.  This  must  now  be  examined  to  evaluate  the  effectiveness  of  the 


shallow  junction  structure.  The  next  chapter  describes  the  analysis. 
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VII  Analysis  of  Results 


This  chapter  presents  the  experimental  results  of  the  preceeding 
chapter  as  curves  and  analyzes  the  performance  of  each  diode  with  varying 
structures.  The  forward  and  reverse  bias  curves  are  displayed  in  Figures 
35  through  39  and  discussed  in  terms  of  diode  characteristics.  The  peak 
optical  characteristics  at  10V  bias  are  plotted  in  Figures  40  through  44 
for  each  set  of  calibrated  data  listed  in  Table  IV.  The  performance  of 
each  device  can  be  compared  to  other  Ge  photodiodes  listed  in  Table  I. 

A  comparison  of  the  I-V  curves  of  the  implanted  and  diffused  diodes 
reveals  the  nicer  overall  characteristics  of  the  diffused  samples.  The 
lower,  flatter  reverse  bias  currents,  characteristic  of  the  diffused 
samples  can  probably  be  attributed  to  the  following  features:  uniformity 
of  the  diffused  junction,  cleaner  processing,  and  mesa  etching.  These 
diodes  were  fabricated  using  well  developed  IC  technology,  i.e.  standard 
masking  and  photoresist  film  processes.  Cleaner  processing  refers  to  the 
fabrication  of  these  devices  in  a  clean  room  using  well-established 
cleaning  procedures.  Due  to  the  very  clean  environment  during  fabrication 
and  tried  masking  procedures,  the  probability  of  a  high  density  of  defects 
is  low,  assuming  a  high  quality  wafer  is  used  to  begin  with.  In  contrast, 
the  implanted  diodes  exhibit  higher  reverse  bias  currents,  3-13  (iA  for 
implanted  diodes  versus  1.5-5 /iA  for  diffused  diodes.  Absorption  problems 
due  to  humidity  are  notoriously  associated  with  Ge  photodiodes.  These 
diodes  were  no  exception.  Under  vacuum  the  higher  (~10mA)  leakage  cur¬ 
rents  of  the  diffused  diodes  was  reduced  to  those  plotted  in  the  Figures. 
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The  implanted  diodes,  on  the  other  hand,  were  not  plagued  by  absorption 
problems.  Implanted  I-V  characteristics  did  not  change  under  vacuum. 

The  purpose  for  mesa-etching  the  two  diffused  and  one  implanted 
diodes  was  to  reduce  the  leakage  current.  The  mesa  structure  was  necessary 
to  define  the  area  of  the  photodiodes  which  had  been  diffused  across  the 
surface  of  the  starting  wafer.  The  decrease  in  leakage  currents  is  sup¬ 
ported  by  comparing  the  I-V  curves  of  the  five  samples.  The  bottom  line 
of  the  low  leakage  current  subject  is  reduced  shot  noise  caused  by  injected 
carriers  in  the  p-n  junction  as  discussed  in  Chapter  II.  A  study  of  noise 
characteristics  was  not  done  in  conjunction  with  this  thesis.  With  a 
"picture"  of  the  five  diodes  in  mind,  an  analysis  of  the  optical  charac¬ 
teristics  of  the  photodiodes  follows. 

The  photoresponse  characteristics  of  each  sample  are  plotted  at  10V 
bias  for  each  set  of  calibration  data,  thereby  making  two  sets  of  optical 
performance  data  available  for  comparison.  In  general,  the  shape  of  the 
response  curves  agrees  with  the  well  known  response  curves  of  Ge  photo¬ 
diodes  as  shown  in  Figure  1  of  Chapter  I.  Peak  response  occurs  at  about 
1.46  /im  for  all  the  samples.  The  drop  in  responsivity  in  the  1.0-1. 3  Mm 
range  is  probably  due  to  the  drop  in  optical  power  at  these  wavelengths 
due  to  reflection  losses.  A  suitable  antireflection  coating  is  likely 
to  solve  this  problem.  A  study  of  the  photoresponse  characteristics  of 
samples  G-7-4  #23 ,  D-6-1  #1,  and  D-6-3  it 6  in  the  Tables  referred  to  in 
Chapter  VI  reveal  responsivities  greater  than  l^tA//UW  and  quantum  effic¬ 
iencies  greater  than  100%  at  peak  performance.  The  problem  may  be  attrib¬ 
uted  to  the  active  areas  of  the  diodes  listed  in  Table  VI  which  are  very 
small  in  comparison  to  the  other  diodes.  The  responsivity  and  quantum 
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efficiency  were  normalized  accordingly  to  compare  the  shapes  of  the  curves, 
which  agree  with  well  known  results. 

Diode  G-7-1  // 9  has  the  best  photocharacteristics  of  the  five  samples. 
It  is  a  simple  planar  implanted  Ge  photodiode  displaying  peak  responsivity 
of  about  0.9/iA//uW  and  peak  efficiency  of  about  80%.  These  performance 
values  can  compete  easily  with  the  other  Ge  photodiodes  described  in 
Table  I.  The  gain  in  performance  is  attributed  to  the  shallow  junction. 

Another  characteristic  of  the  Ge  photodiodes  is  carrier  lifetime 
which  plays  a  role  in  the  denominator  of  the  leakage  current  as  shown  in 
Equation  (26)  (Ref  1;  4;  7;  18;  19;  30;  33;  43;  47;  55).  An  increase  in 
carrier  lifetime  would  decrease  the  leakage  current.  In  contrast  to 
Ando's  planar  avalanche  photodiode  in  Table  I  which  claims  a  lifetime  of 
0.9ns,  the  measured  lifetime  for  the  diffused  diodes  is  2.25 /is;  there 
is  an  obvious  discrepancy  in  the  orders  of  magnitudes.  High  responsivity 
in  the  APD  is  due  to  the  photomultiplication  mechanism  and  requires  low 
values  of  lifetime  to  increase  the  speed  performance.  On  the  other  hand, 
in  the  p-n  junction,  the  high  response  is  due  to  the  shallow  junction 
structure  and  the  "instantaneous"  collection  of  photogenerated  carriers 
in  a  region  where  there  is  a  high  amount  of  light  absorbed.  The  longer 
lifetime  reduces  the  leakage  current  responsible  for  noise  in  p-n  diodes, 
whereas  the  low  lifetime  of  the  APD's  makes  them  inherently  noisy  devices. 

In  summary,  the  response  of  the  planar  Ge  photodiodes  displayed  peak 
performance  at  the  expected  1.46/zm  wavelength  with  competitive  values 
of  photoresponse  and  quantum  efficiency.  The  use  of  an  AR  coating  is 
also  expected  to  improve  the  photocharacteristics.  Better  IC  fabrication 
facilities  would  be  expected  to  improve  the  performance  even  more. 
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Fig  35.  I-V  characteristics  of  G-7-1//9. 
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Fig  37. 


T'V  characterlstiCs 
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VOLTAGE  (V) 


1-V  characteristics  of  D-6-1//6. 


1.0 


Fig  40(a).  Photoresponse  and  quantum  efficiency  curves  of  G— 7— 1//9 
at  10  V  bias  for  H*  calibration  data. 
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WAVELENGTH  (tfm) 

Fig  41(a).  Photoresponse  and  quantum  efficiency  curves  of  G-7-4//9 
at  10  V  bias  for  H*  calibration  data. 
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Fig  43(b).  Photoresponse  and  quantum  efficiency  curves  of  D' 
at  10  V  bias  for  H**  calibration  data. 
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VIII  Conclusions  and  Recommendations 

In  conclusion,  the  feature  which  increases  the  performance  character¬ 
istics  of  the  everyday  Ge  photodiode  is  the  implementaion  of  a  shallow 
junction  to  "capture"  most  of  the  light  absorbed  within  5 fim  from  the 
surface.  The  best  response  values  for  an  impalnted  sample  were  R=0.9jzA/MW 
and  7)  =  80%,  competitive  with  other  diffused  Ge  photodiode  structures 
which  have  been  studied.  The  performance  can  be  further  improved  with  an 
AR  coating  to  reduce  leakage  currents,  and  encapsulant  to  avoid  problems 
with  humidity,  guard  rings,  side  contacts,  and  other  conventional  improve¬ 
ments  not  feasible  in  research  facilities. 

In  the  area  of  Ge-Si  alloy  crystal  growth  which  is  in  its  very  early 
stages  of  development,  a  quantum  mechanical  analysis  of  the  Ge-Si  energy 
band  structure  should  be  made  using  the  tight-binding  approximation.  A 
repeated  application  of  the  LCAO  method  as  outlined  by  Slater  and  Koster 
should  be  applied  to  a  representative  group  of  unit  cells  whose  total 
volume  represents  a  concentration  10%  Si  and  90%  Ge.  The  band  structure 
of  the  Ge-Si  alloy  would  be  a  weighted  average  of  the  energy  band  structures 
of  the  individual  unit  cells.  The  results  of  this  modeling  should  pre¬ 
dict  the  energy  band  structure  of  the  crystal  alloy  and  in  doing  so  its 
effective  mass  and  energy  gap,  thereby  predicting  the  behavior  of  electrons 
and  holes  under  the  influence  of  an  applied  field.  Also  recommended  are 
optical  and  cyclotron  resonance  experiments  to  quantify  the  theoretical 
results  for  bandgap  and  effective  mass  values  respectively.  These  values, 
at  crucial  points  in  the  symmetry  points  of  the  Brillouin  zone,  would  be 
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helpful  in  fitting  the  constants  arising  in  the  LCAO  method  at  nonsymmetry 
points  in  the  zone.  Once  the  Ge-Si  alloy  is  grown  and  is  available  in 
wafer  form,  a  shallow-junction  Ge-Si  photodiode  should  be  fabricated  and 
tested.  It  is  expected  that  Ge-Si  photodetector  will  exhibit  a  peak 
photoresponse  of  about  0.8  -  0.9/^A/yUW  at  1.27/im  where  the  optical 
fiber  experiences  minimum  attenuation  and  dispersion,  thereby  optimizing 
system  performance  in  fiber-optic  applications. 
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APPENDIX  A 


Synchronous  Detection 

A  Lock-In  Voltmeter  is  substituted  for  the  amplifier  in  Figure  30  and 
a  chopper  blade  is  inserted  between  the  source  and  the  photodetector 
(Ref  6;  26;  41;  50).  The  chopper  blade  rotates  at  constant  speed  and  the 
detector  produces  a  fluctuating  voltage^^at  the  input  to  the  Lock-In 
Voltmeter  as  shown  in  Figure  45.  In  the  Lock-In  Voltmeter,  the  DC  part 
of  the  signal  is  rejected.  The  remaining  noise  and  square  components  are 
amplified  which  causes  high  level  push-pull  signals  to  be  generated  by  the 
transformer , (C^and (5^.  The  chopper  blade  is  operated  synchronously  with 
the  push-pull  signals.  The  positive  half  cycles  of  the  push-pull  signals 
are  mixed  producing  a  full-wave  rectified  output  (e)  upon  which  a  the 
noise  is  superimposed.  The  random  noise  almost  cancels  and  an  RC  filter 
smooths  out  the  DC  output,  the  output  of  the  Lock-In  Voltmeter,  (?)  . 

The  synchronous  detector  system  is  almost  free  from  zero  error  and  greater 
accuracy  of  measurements  is  obtained. 
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Fig  45.  Application  and  operation  of  the  lock- in-amplifier  in  a  synchronous 
detection  scheme  (Ref  6;  26). 


APPENDIX  B 


Experimental  Procedure 
Mounting  the  Chip  for  Testing 

The  steps  for  mounting  a  chip  for  testing  can  be  followed  using 
Figure  33  and  arc  described  here  in  more  detail.  Leads  were  soldered  to 
the  appropriate  pins:  6,  9,  20,  23,  and  15.  The  chip  was  placed  vertic¬ 
ally  on  the  vacuum  header  and  clamped  facing  the  output  slit  of  the  mono¬ 
chromator  as  shown  in  Figure  46,  i.e.  in  direct  alignment  with  the  light 
path  of  the  radiant  source.  The  cylindrical  vacuum  chamber  was  slipped 
over  the  header  aligning  the  silicon  glass  input  window  with  the  window 
of  the  monochromator  so  that  the  sample  could  receive  incident  light. 

The  leads  from  the  chip  were  pulled  through  the  circular  hole  in  the  back 
of  the  vacuum  chamber  as  in  Figure  47.  Each  labeled  lead  was  inserted 
into  its  matching  hole  in  the  5-pin  socket  of  the  connector  with  external 
coaxial  output  pins  as  shown  in  Figure  48.  The  connector  was  held  in 
place  with  four  screws  and  the  "0"  ring  on  the  vacuum  chamber  formed  the 
seal  for  the  vacuum.  The  air-intake  valve  at  the  base  of  the  vacuum 
chamber  was  securely  closed  and  the  main  valve  to  the  vacuum  pump  was 
opened  to  allow  a  vacuum  to  build  up  inside  the  chamber.  The  vacuum  al¬ 
lowed  the  diode  to  operate  in  a  protective  environment  without  the  influ¬ 
ence  of  degrading  factors  in  the  ambient,  such  as  dust  and  humidity. 

Recall  that  these  diodes  were  not  sealed  in  any  way.  Vacuum  was  allowed 
-4 

to  reach  10  torr  (0.1  ^m  Hg) . 

Once  the  sample  was  mounted  in  the  header,  a  vacuum  developed  ,  and 
the  globar  warmed  up  to  operating  specifications  (230  watts)  photoresponse 


back  of  the  header.  The  diode  faces  the 
window  of  the  monochromator. 


Fig  48.  Photo  showing  the  leads  from  the  diode  being  inserted  into  a 
5-pin  socket  of  the  connector  with  coaxial  output  pins. 


measurements  and  diode  characteristics  could  be  taken  on  the  device. 
Coaxial  cables  were  used  to  connect  each  diode  to  a  curve  tracer  to  obtain 
the  forward  and  reverse  bias  diode  characteristics.  A  black  cloth  was 
draped  over  the  header  of  the  vacuum  chamber  to  prevent  ambient  light  from 
filtering  into  the  chamber  and  to  prevent  the  source  from  illuminating 
the  sample  during  the  check  for  device  characteristics.  The  black  cloth 
and  coaxial  cables  were  then  removed  from  the  header  to  perform  photo¬ 
response  measurements. 

The  cables  from  the  preamplifier  were  attached  to  the  connector  to 
engage  one  of  the  photodiodes  into  the  system(e.g.  #23  and  #15). 
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By  turning  the  drum  of  the  monochromator,  it  was  possible  to  manually  drive 
the  optical  prisms  of  the  system  so  that  the  output  radiation  of  the  glo- 
bar  emitted  light  through  a  spectrum  of  wavelengths.  These  ranged  from 
1.0/im  to  2.03  ^im.  The  silicon  window  on  the  vacuum  header  started  to 
cut  off  wavelengths  below  l.O^Lim,  that  is,  below  l.Oyum  there  was  not 
100%  transmission. 

Equipment  Turn-On  and  Test  Sequence 

Before  doing  anything,  the  cooling  water  for  the  globar  was  turned  on. 
The  temperature  dial  at  the  top  of  the  rack  was  set  to  room  temperature, 
290°K.  The  system  has  cooling  capabilities,  but  these  Ge  diodes  were 
meant  to  operate  at  room  temperature.  The  globar  power  was  checked  to 
make  sure  the  dial  was  turned  all  the  way  to  the  left,  to  zero.  The  power 
switch  to  the  globar  was  turned  on,  set  to  20V,  and  allowed  to  warm  up 
slowly  in  50-75  watt  increments  to  230  watts,  where  the  source  spectral 
characteristics  are  calibrated. 

The  oscillator  and  lock-in  voltmeter  were  turned  on.  The  KEPCO 
bias  supply  was  set  to  0  volts  before  being  turned  on.  The  digital  volt¬ 
meter  was  used  as  a  check  to  monitor  the  bias  voltage.  The  digital  record¬ 
er  was  set  to  standby  and  turned  on.  The  function  of  the  programmer  was 
to  set  the  speed  of  the  chopper  to  139  Hz.  This  frequency  was  chosen  as 
the  frequency  with  the  lowest  susceptibility  to  noise  pick-up.  The 
programmer  wal  also  set  for  the  number  of  openings  on  the  chopper  wheel, 
six  in  this  case. 

The  x-y  staging  was  adjusted  for  maximum  output.  The  maximum  peak- 

to-peak  output  was  observed  on  the  oscilloscope  and  is  explained  using 

Figure  49.  To  synchronize  the  chopper  speed  with  the  lock-in  voltmeter, 
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the  latter  was  set  to  the  calibrate  mode.  The  chopper  speed  was  trimmed 
to  align  the  scope  display  to  get  a  picture  of  a  full-rectified  wave  as 
in  Figure  49.  Since  the  output  of  the  Brower  is  always  0-10  volts,  cal¬ 
ibration  was  required  for  each  scale.  For  example,  if  the  input  signal 
was  between  0  and  25  mV,  the  maximum  reading  on  the  voltmeter  would  be 
4  volts.  More  simply,  the  scaling  is: 

0  <  input  signal  <  25  mV 

0  <  Brower  output  <  10  V 

input  signal  =  25  mV  x  (Brower  output) 

10  V 

The  gain  on  the  lock-in  voltmeter  was  turned  to  adjust  the  Brower  output 
to  its  maximum  and  the  calibration  switch  was  turned  to  normal.  To  find 
the  maximum  p-p  output,  the  phasing  on  the  lock-in  voltmeter  is  similar 
to  the  adjustment  shown  in  Figure  49. 

To  take  data,  the  programmer  was  put  on  hold  and  the  data  recorder 
on  operate.  Moving  from  high  to  low  drum  numbers  in  the  back  of  the 
monochromator  (so  that  the  optics  would  work  properly)  as  shown  in 
Figure  50,  the  gate  of  the  programmer  was  triggered  by  pushing  the  re¬ 
mote  button  to  take  data  when  ready.  At  the  end  of  a  run,  the  5055 
recorder  was  put  on  standby  and  the  3490  programmer  was  taken  off  hold 
to  sample  voltages  at  regular  intervals. 

Removing  the  Sample  From  the  Evacuated  Chamber. 

After  completion  of  the  tests  for  a  particular  chip,  the  valve  on 
the  main  line  to  the  vacuum  pump  was  closed  to  isolate  the  pump  from  the 
chamber.  Air  was  admitted  into  the  chamber  by  a  valve  at  its  base. 
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BEFORE 


AFTER 


Fig  50.  The  rotating  drum  and  remote  switch  used  to  scan  the  spectral 
output  of  the  globar  and  trigger  the  data  recorder. 


The  screws  were  removed  from  the  measuring  head  which  was  carefully  re¬ 
moved  after  separating  the  leads  from  the  5-pin  socket.  The  vacuum 
chamber  was  removed  by  pulling  straight  up  and  the  sample  was  removed 
from  the  holder  by  unscrewing  the  clamping  screw. 

Equipment  Turn-Off  Procedure. 

All  bias  was  removed  from  the  diode  after  a  test  was  completed;  the 
KEPCO  bias  supply  was  turned  off  and  set  to  zero  as  a  precaution.  The 
lock-in  voltmeter,  oscillator,  and  chopper  wheel  were  turned  off  as  well 
as  all  voltmeters  and  the  digital  recorder.  The  globar  power  was  turned 
down  slowly  in  50  watt  increments  to  prolong  the  life  of  the  globar. 

Once  the  globar  power  was  zero,  it  was  switched  off,  and  the  system 
power  turned  off.  The  final  step  was  to  close  the  water  valve  for 
cooling  the  globar. 
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APPENDIX  C 


I-V  Points  Used  to  Plot  the  Diode  Characteristics  of  Each  Sample 


TABLE  C-l 


Diode 

Characteristics 

of  G-7-1  #9 

Reverse  Bias 

Forward 

Bias 

V(V) 

I(UA) 

l - - - - - 

V(V) 

I(UA) 

V(V) 

igwA) 

0.05 

3.3 

2 

i 

6.5 

0.05 

7.5 

0.10 

4.5 

3 

7.0 

0.10 

17.0 

0.15 

4.8 

4 

7.5 

0.15 

30.0 

0.20 

5.0 

5 

8.0 

0.20 

40.0 

0.25 

5.1 

10 

10.0 

1.0 

250.0 

0.30 

5.2 

15 

12.0 

0.4 

5.3 

20 

14.0 

0.5 

5.4 

25 

16.0 

1.0 

5.8 

A 
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Forward  Bias 


Kj/A) 

V(V) 

I(/iA) 

17.5 

0.05 

18.0 

18.0 

0.06 

23.0 

18.0 

0.07 

28.0 

19.0 

0.08 

36.0 

20.0 

0.09 

46.0 

21.0 

0.10 

52.5 

21.5 

0.15 

95.0 

22.0 

0.20 

150.0 

22.5 

0.25 

225.0 

24.0 

0.30 

320.0 

28.5 

0.35 

400.0 

r 
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Reverse  Bias 


Forward  Bias 


V(V) 

Ify/A) 

V(V) 

I0*A) 

V(V) 

IOUA) 

0.05 

2.75 

27.5 

6.0 

0.01 

1.0 

0.10 

3.1 

30.0 

6.3 

0.02 

2.0 

0.15 

4.6 

32.5 

7.0 

0.03 

3.6 

0.20 

5.0 

35.0 

8.0 

0.04 

5.6 

0.25 

5.2 

37.5 

10.0 

0.05 

7.5 

0.30 

5.3 

40.0 

20.0 

0.06 

9.5 

0.35 

5.35 

0.07 

13.0 

0.40 

5.4 

0.08 

18.0 

0.45 

5.4 

0.09 

24.0 

4.5 

5.5 

0.10 

33.0 

20.0 

5.6 

0.15 

135.0 

25.0 


5.8 


0.25 


660.0 


APPENDIX  D 


Output  Voltages  Recorded  Across  the  LOOK  Resistor  for  Each  Test  Run 


TABLE  D-l 


Output  Photovoltage  of  G-7-I  #9  for  H*  Calibration  Data 


X 

Cum) 

0V  Bias 

V(//V) 

5V  Bias 

V(UV) 

10V  Bias 
V(^V) 

1.04 

170 

2625 

2775 

1.11 

1256 

19450 

19950 

1.19 

3057 

48950 

51125 

1.28 

4517 

74300 

76175 

1.37 

5794 

97150 

98425 

1.46 

7173 

123550 

125075 

1.55 

6965 

119925 

122300 

1.65 

5388 

90100 

91325 

1.75 

4214 

68400 

69575 

1.85 

2567 

41475 

41775 

2.03 

114 

1950 

1925 

TABLE  D-2 

Output  Photovoltage  of  G-7-1  #9 

for  H**  Calibration 

Data 

X 

(Um) 

0V  Bias 

V( fjV) 

5V  Bias 

VC UV) 

10V  Bias 
V(f/V) 

1.04 

188 

2825 

2950 

1.11 

1282 

20300 

20700 

1.19 

2976 

50350 

50125 

1.28 

4296 

75950 

75025 

1.37 

5452 

98900 

97675 

1.46 

6596 

122225 

120550 

1.55 

6342 

117150 

115725 

1.65 

6342 

86800 

86850 

1.75 

3948 

66450 

65950 

1.85 

2622 

42600 

42000 

2.03 

134 

2000 

2150 

( 


134 


TABLE  D-3 


Output  Photovoltage  of  G-7-4  #23 

for  H*  Calibration 

Data 

X 

0V  Bias 

5V  Bias 

10V  Bias 

(fjm) 

V(MV) 

V(/iV) 

V(/iV) 

1.04 

1185 

8350 

8175 

1.11 

6988 

54275 

54000 

1.19 

14343 

135275 

133450 

1.28 

19570 

202250 

199925 

1.37 

23718 

260600 

260500 

1.46 

27610 

327750 

324700 

1.55 

25853 

295600 

298450 

1.65 

19383 

192525 

197025 

1.75 

14300 

126700 

125450 

1.85 

8063 

66625 

67300 

2.03 

318 

2725 

2775 

TABLE  D-4 

Output  Photovoltage  of  G-7-4  #23 

for  H**  Calibration 

Data 

X 

(fjm) 

0V  Bias 

VQ/V) 

5V  Bias 

V(//V) 

10V  Bias 
VCuV) 

1.04 

655 

11050 

11600 

1.11 

4065 

72500 

76100 

1.19 

8950 

174500 

179800 

1.28 

12465 

261700 

270150 

1.37 

15395 

338450 

349600 

1.46 

17948 

419250 

431100 

1.55 

16750 

379000 

387650 

1.65 

11978 

246450 

250800 

1.75 

8273 

159450 

163650 

1.85 

4575 

81050 

82900 

2.03 

168 

2750 

2800 
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TABLE  D-5 


Output  Photovoltage  of  G-8-1  if 23  for  H*  Calibration  Data 


X 

</im) 

0V  Bias 

V(flV ) 

5V  Bias 

V(f/V) 

10V  Bias 
V(/iV) 

1.04 

216 

5300 

5575 

1.11 

1441 

36450 

38925 

1.19 

3573 

95200 

100075 

1.28 

5500 

148725 

157850 

1.37 

7329 

203975 

212400 

1.46 

9459 

256525 

281050 

1.55 

9838 

281525 

296500 

1.65 

8135 

228450 

239200 

1.75 

5846 

162825 

171825 

1.85 

3531 

94800 

99125 

2.03 

195 

4975 

5025 

TABLE  D-6 

Output  Photovoltage  of  G-8-1  if 23 

for  H**  Calibration 

Data 

X 

0V  Bias 

5V  Bias 

10V  Bias 

(flm) 

V(/UV) 

V(/iV) 

V(/7V' 

1.04 

165 

3550 

4010 

1.11 

1261 

30600 

31900 

1.19 

3205 

81250 

1.28 

4995 

130950 

137400 

1.37 

6603 

178400 

184550 

1.46 

8396 

234100 

244650 

1.55 

8838 

249250 

258550 

1.65 

7459 

203050 

211700 

1.75 

5544 

145900 

148700 

1.85 

3171 

80850 

83550 

2.03 

156 

3550 

3750 

136 


TABLE  D-7 


Output  Photovoltage  of  D-6-1  #6  for  H*  Calibration  Data 


X 

(pm) 

0V  Bias 

V(//V) 

5V  Bias 

V(jlV) 

10V  Bias 
V(//V) 

1.04 

830 

3525 

1425 

1.11 

5565 

26975 

22375 

1.19 

12995 

71025 

59125 

1.28 

18515 

110850 

91550 

1.37 

23055 

148625 

122100 

1.46 

28080 

195950 

162500 

1.55 

28205 

201025 

167650 

1.65 

22895 

152050 

126975 

1.75 

16650 

102950 

84500 

1.85 

8895 

48400 

41550 

2.03 

360 

1650 

1725 

TABLE  D-8 

Output  Photovoltage  of  D-6-1  # 6 

for  H**  Calibration 

Data 

X 

(Jim) 

0V  Bias 

V(J1V) 

5V  Bias 

V(flV) 

10V  Bias 
V(j«V) 

1.04 

755 

5000 

4450 

1.11 

4695 

35250 

32700 

1.19 

10470 

87800 

81450 

1.28 

14660 

131900 

723100 

1.37 

18080 

173650 

161500 

1 .46 

21375 

218950 

203750 

1.55 

21145 

221450 

206100 

1.65 

17010 

170150 

160150 

1.75 

12370 

116750 

108200 

1.85 

7260 

62950 

58500 

2.03 

345 

2600 

2050 
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TABLE  D-9 


Output  Photovoltage  of  D-6-3  #6  for  H*  Calibration  Data 


A 

(jim) 

0V  Bias 

V(/1V) 

5V  Bias 

V(J/V) 

10V  Bias 
V0/V) 

1.04 

1090 

10000 

10350 

1.11 

6580 

64500 

64650 

1.19 

14407 

152200 

154400 

1.28 

20235 

227300 

229950 

1.37 

24790 

294800 

297500 

1.46 

29217 

370450 

374050 

1.55 

28597 

365950 

368800 

1.65 

22277 

279000 

280750 

1.75 

15997 

192800 

194850 

1.85 

8992 

99800 

102050 

400 

4200 

4250 

TABLE  D-10 

Output  Photovoltage  of  D-6-3  #6 

for  H**  Calibration 

Data 

A 

(flm) 

0V  Bias 

V(pCV) 

5V  Bias 

V(fiV) 

10V  Bias 
V(/iV) 

1.04 

470 

10650 

11050 

1.11 

3203 

73350 

72700 

1.19 

7328 

176000 

175400 

1.28 

10673 

263250 

261400 

1.37 

13340 

343400 

340400 

1.46 

15908 

424200 

414300 

1.55 

15238 

414500 

422850 

1.65 

11125 

310900 

309200 

1.75 

7573 

214100 

214350 

1.85 

4243 

116550 

115750 

2.03 

158 

4450 

4500 
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APPENDIX  E 


Photocurrent,  Incident  Optical  Power,  Responsivity,  and  Quantum 
Efficiency  of  Each  Sample  per  Test  Run 


TABLE  E-l 

Photocharacteristics  of  G-7-1  #9  at  O.lmV  for  H*  Calibration  Data 


X 

I 

H 

R 

V 

<//m) 

(Uk) 

(jUW) 

(tik/flM) 

(%) 

1.04 

0.00170 

0.059 

0.029 

3.4 

1.11 

0.01256 

0.295 

0.043 

4.8 

1.19 

0.03057 

0.648 

0.047 

4.9 

1.28 

0.04517 

0.912 

0.050 

4.8 

1.37 

0.05794 

1.104 

0.052 

4.8 

1.46 

0.07173 

1.357 

0.053 

4.5 

1.55 

0.06965 

1.521 

0.046 

3.7 

1.65 

0.05388 

1.616 

0.033 

2.5 

1.75 

0.04214 

1.640 

0.026 

1.8 

1.85 

0.02567 

1.583 

0.162 

1.1 

2.03 

0.00114 

1.621 

0.001 

TABLE  E-2 

Photocharacteristics 

of  G-7-1  #9  at 

5V  for  H*  Calibration 

Data 

X 

I 

H 

R 

V 

(A/m) 

(//A) 

<A/w) 

(Uk/fj 

(%) 

1.04 

0.02625 

0.059 

0.445 

53.0 

1.11 

0.19450 

0.295 

0.659 

73.7 

1.19 

0.48950 

0.648 

0.755 

78.7 

1.28 

0.74300 

0.912 

0.815 

78.9 

1.37 

0.97150 

1.104 

0.880 

79.6 

1.46 

1.23550 

1.357 

0.910 

77.3 

1.55 

1.19925 

1.521 

0.788 

63.1 

1.65 

0.90100 

1.616 

0.558 

41.9 

1.75 

0.68400 

1.640 

0.417 

29.6 

1.85 

0.41475 

1.583 

0.262 

17.6 

2.03 

0.01950 

1.621 

0.120 

7.3 

i 

4 


TABLE  E-3 


Photocharacteris : ..  ,  of  G-7-1  #9  at  10V  for  H*  Calibration  Data 


X 

I 

H 

R 

V 

(jm) 

VIA) 

V/w) 

(/!A//iW) 

(%) 

1.04 

0.02775 

0.059 

0.470 

56.1 

1.11 

0.19950 

0.295 

0.676 

75.5 

1.19 

0.51125 

0.648 

0.789 

82.2 

1.28 

0.76175 

0.912 

0.835 

80.9 

1.37 

0.98425 

1.104 

0.892 

80.7 

1.46 

1.25075 

1.357 

0.922 

78.3 

1.55 

1.22300 

1.521 

0.804 

64.3 

1.65 

0.91325 

1.616 

0.565 

42.5 

1.75 

0.69575 

1.640 

0.424 

30.1 

1.85 

0.41775 

1.583 

0.264 

17.7 

2.03 

0.01925 

1.621 

0.012 

0.7 

TABLE  E-4 

Photocharacteristics 

of  G-7-1  //9  at 

0V  for  H**  Calibration 

Data 

X 

I 

li 

R 

V 

(fim) 

(MA) 

(/1W) 

QlA/jiW) 

(%) 

1.04 

0.00188 

0.0464 

0.041 

4.8 

1.11 

0.01282 

0,2998 

0.043 

4.8 

1.19 

0.02976 

0.6460 

0.046 

4.8 

1.28 

0.04296 

0.9225 

0.047 

4.5 

1.37 

0.05452 

1.0893 

0.050 

4.5 

1.46 

0.06596 

1.3152 

0.050 

4.3 

1.55 

0.06342 

1.4671 

0.043 

3.5 

1.65 

0.04886 

1.5410 

0.032 

2.4 

1.75 

0.03948 

1,5664 

0.025 

1.8 

1.85 

0.02622 

1.5057 

0.017 

1.2 

2.03 

0.00134 

1.5368 

0.001 
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TABLE  E-5 


Photocharacteristics  of  G-7-1  // 9  at  5V  for  H**  Calibration  Data 


X 

I 

H 

R  . 

»7 

(/An) 

(/UA) 

<JM) 

(/M//iW) 

(%) 

1.04 

0.02825 

0.0464 

0.615 

72.5 

1.11 

0.20300 

0.2998 

0.677 

75.6 

1.19 

0.50350 

0.6460 

0.779 

81.2 

1.28 

0.75950 

0.9225 

0.823 

79.7 

1.37 

0.98900 

1.0893 

0.908 

82.2 

1.46 

1.22225 

1.3152 

0.929 

78.9 

1.55 

1.17150 

1.4671 

0.799 

63.9 

1.65 

0.86800 

1.5410 

0.563 

42.3 

1.75 

0.66450 

1.5664 

0.424 

30.0 

1.85 

0.42600 

1.5051 

0.283 

19.0 

2.03 

0.02000 

1.5368 

0.013 

.79 

TABLE  E-6 

Photocharacteristics  of  G-7-1  1f9  at  10V  for  H**  Calibration  Data 


X 

(/An) 

I 

(//A) 

H 

<//w) 

R 

(UA///W) 

V 

(%) 

1.04 

0.02950 

0.0464 

0.636 

75.8 

1.11 

0.20700 

0.2998 

0.695 

77.1 

1.19 

0.50125 

0.6460 

0.776 

80.9 

1.28 

0.75025 

0.9225 

0.813 

78.8 

1.37 

0.97675 

1.0893 

0.897 

81.2 

1.46 

1.20550 

1.3152 

0.917 

77.8 

1.55 

1.15725 

1.4671 

0.789 

63.1 

1.65 

0.86850 

1.5410 

0.564 

42.4 

1.75 

0.65950 

1.5664 

0.438 

31.0 

1.85 

0.42000 

1.5051 

0.279 

18.7 

2.03 

0.02150 

1.5368 

0.014 

0.9 

TABLE  E-7 


Photocharacteristics  of  G-7-4  // 23  at  OV  for  H*  Calibration  Data 


X 

I 

H 

R 

V 

CM  in) 

(MA) 

(MW) 

(MA/mw) 

(%) 

1.04 

0.00655 

0.0285 

0.230 

27.4 

1.11 

0.04065 

0.1840 

0.221 

24.7 

1.19 

0.08950 

0.3966 

0.226 

23.5 

1.28 

0.12465 

0.5664 

0.220 

21.3 

1.37 

0.15395 

0.6687 

0.230 

20.8 

1.46 

0.17948 

0.8074 

0.222 

18.9 

1.55 

0.16750 

0.9007 

0.186 

14.9 

1.65 

0.11978 

0.9461 

0.127 

9.5 

1.75 

0.08273 

0.9616 

0.086 

6.1 

1.85 

0.04575 

0.9241 

0.050 

3.4 

2.03 

0.00168 

0.9435 

0.002 

.12 

TABLE  E-8 

Photocharacteristics  of 

G-7-4  //23  at 

5V  for  H* 

Calibration 

Data 

X 

I 

H 

R 

R' 

V 

V' 

Oim) 

(MA) 

(MW) 

OiA/MW) 

(MA/mw) 

(%) 

(%) 

1.04 

0.11050 

0.0285 

0.3877 

0.7466 

46.2 

100.0 

1.11 

0.72500 

0.1840 

0.3940 

0.7587 

44.0 

95.2 

1.19 

1.74500 

0.3966 

0.4400 

0.8473 

45.8 

99.1 

1.28 

2.61700 

0.5664 

0.4620 

0.8897 

44.8 

97.0 

1.37 

3.38450 

0.6687 

0.5061 

0.9746 

45.8 

99.1 

1.46 

4.19250 

0.8074 

0.5193 

1.0000 

44.1 

96.1 

1.55 

3.79000 

0.9007 

0.4694 

0.9039 

39.9 

86.4 

1.65 

2.46450 

0.9461 

0.2605 

0.5016 

19.6 

42.4 

1.75 

1.59450 

0.9616 

0.1658 

0.3193 

11.7 

25.3 

1.85 

0.81050 

0.9241 

0.0877 

0.1689 

5.9 

12.8 

2.03 

0.02750 

0.9435 

0.029 

0.0558 

0.2 

0.4 

TABLE  E-9 


Photocharacteristics  of  G-7-4  #23  at  10V  for  H*  Calibration  Data 


X 

I 

H 

R 

R' 

V 

r 

(m®) 

(ma) 

</iW) 

(/iA/fiW) 

(ma/^w) 

<%) 

(%> 

1.04 

0.11600 

0.0285 

0.4070 

0.7623 

48.5 

100.0 

1.11 

0.76100 

0.1840 

0.4136 

0.7747 

46.2 

95.3 

1.19 

1.79800 

0.3966 

0.4533 

0.8490 

hi. 2 

97.3 

1.28 

2.70150 

0.5664 

0.4770 

0.8934 

46.2 

95.3 

1.37 

3.49600 

0.6687 

0.5228 

0.9792 

47.3 

97.5 

1.46 

5.31100 

0.8074 

0.5339 

1.0000 

45.3 

93.4 

1.55 

3.87650 

0.9007 

0.4304 

0.8061 

34.4 

70.9 

1.65 

2.50800 

0.9461 

0.2651 

0.4965 

19.9 

41.0 

1.75 

1.63650 

0.9616 

0.1702 

0.3188 

12.1 

24.9 

1.85 

0.82900 

0.9241 

0.0897 

0.1680 

6.0 

*  12.4 

2.03 

0.02800 

0.9435 

0.0030 

0.0056 

.2 

0.4 

TABLE  E-10 

Photocharacteristics  of  G-7-4  #23  at  O.lmV  for  H**  Calibration  Data 


X 

I 

H 

R 

V 

(n®) 

(ma) 

(MW) 

(jiA/nW) 

(%) 

1.04 

0.01185 

0.0356 

0.332 

39.7 

1.11 

0.06988 

0.1778 

0.393 

43.9 

1.19 

0.14343 

0.3899 

0.368 

38.3 

1.28 

0.19570 

0.5486 

0.357 

34.6 

1.37 

0.23718 

0.6642 

0.357 

32.3 

1.46 

0.27610 

0.8166 

0.338 

28.7 

.55 

0.25853 

0.9157 

0.282 

22.6 

1.65 

0.19383 

0.9728 

0.199 

15.0 

1.75 

0.14300 

0.9868 

0.145 

10.3 

1.85 

0.08063 

0.9525 

0.085 

5.7 

2.03 

0.00318 

0.9754 

0.003 

0.2 
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TABLE  E-ll 


Photocharacteristics  of  G-7-4  #23  at  5V  for  H**  Calibration  Data 


X 

I 

H 

R 

R’ 

ri 

V' 

(M) 

(MW) 

(MA/mW) 

(nA/fiW) 

(%) 

(%) 

1.04 

0.08350 

0.0356 

2.345 

0.584 

280 

77 

1.11 

0.54275 

0.1778 

3.053 

0.761 

341 

94 

1.19 

1.35275 

0.3899 

3.469 

0.864 

362 

100 

1.28 

2.02250 

0.5486 

3.687 

0.919 

357 

99 

1.37 

2.60600 

0.6642 

3.924 

0.978 

355 

98 

1.46 

3.27750 

0.8166 

4.014 

1.000 

341  . 

94 

1.55 

2.95600 

0.9157 

3.228 

0.804 

258 

71 

1.65 

1.92525 

0.9728 

1.979 

0.493 

149 

41 

1.75 

1.26700 

0.9868 

1.284 

0.320 

90 

25 

1.85 

0.66625 

0.9525 

0.699 

0.174 

47 

13 

2.03 

0,02725 

0.9754 

0.028 

0.007 

2 

1 

TABLE  E-12 

Photocharacteristics  of  G-7-4  #23  at  10V  for  H**  Calibration  Data 


X 

1 

H 

R 

R' 

V 

T}’ 

OiA) 

(MW) 

(ma/mw) 

(ma/mw) 

(%) 

(%) 

1.04 

0.08175 

0.0356 

2.296 

0.577 

274 

77 

1.11 

0.54000 

0.1778 

3.037 

0.764 

339 

95 

1.19 

1.33450 

0.3899 

3.423 

0.861 

357 

100 

1.28 

1.99925 

0.5486 

3.644 

0.916 

353 

99 

1.37 

2.60500 

0.6642 

3.922 

0.986 

355 

99 

1.46 

3.24700 

0.8166 

3.976 

1.000 

338 

95 

1.55 

2.98450 

0.9157 

3.259 

0.820- 

261 

73 

1.65 

1.97025 

0.9728 

2.025 

0.509 

152 

43 

1.75 

1.25450 

0.9868 

1.271 

0.320 

90 

25 

1.85 

0.67300 

0.9525 

0.707 

0.118 

47 

13 

2.03 

0.02775 

0.9754 

0.028 

0.007 

2 

1 

TABLE  E-13 


Photocharacteristics  of  G-8-1  #23  at  OV  for  H*  Calibration  Data 


X 

(pm) 

I 

( pA ) 

H 

(pw) 

R 

(pA/pW) 

V 

(%) 

1.04 

0.00216 

0.181 

0.012 

1.4 

1.11 

0.01441 

0.907 

0.016 

1.8 

1.19 

0.03573 

1.989 

0.018 

1.9 

1.28 

0.05500 

2.799 

0.020 

1.9 

1.37 

0.0^329 

3.389 

0.022 

2.0 

1.46 

0.09459 

4.167 

0.023 

1.9 

1.55 

0.09838 

4.672 

0.021 

1.7 

1.65 

0.08135 

4.964 

0.016 

1.2 

1.75 

0.05846 

5.035 

0.012 

0.8 

1.85 

0.03531 

4.860 

0.007 

0.5 

2.03 

0.00195 

4.977 

TABLE  E-14 

Photocharacteristics  of  G-8-1  #23  at  5V  for  H*  Calibration  Data 


X 

(pm) 

I 

(PA) 

H 

(pW) 

R 

(pA/pW) 

TJ 

(%) 

1.04 

0.05300 

0.181 

0.293 

34.9 

1.11 

0.36450 

0.907 

0.402 

44.9 

1.19 

0.95200 

1.989 

0.479 

49.9 

1.28 

1.48725 

2.799 

0.531 

51.5 

1.37 

2.03975 

3.389 

0.602 

54.5 

1.46 

2.56525 

4.167 

0.616 

52.3 

1.55 

2.81525 

4.672 

0.603 

48.2 

1.65 

2.28450 

4.964 

0.460 

34.6 

1.75 

1.62825 

5.035 

0.323 

22.9 

1.85 

0.94800 

4.860 

0.195 

13.1 

2.03 

0.04975 

4.977 

0.010 

0.6 

( 
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TABLE  E-15 


Photocharacteristics  of  G-8-1  #23  at  10V  for  H*  Calibration  Data  ! 


X 

I 

H 

R 

7? 

(pm) 

(pA) 

(pW) 

(pA/pW) 

(%) 

1.04 

0.05575 

0.181 

0.308 

36.7 

1.11 

0.38925 

0.907 

0.429 

47.9 

1.19 

1.00075 

1.989 

0.503 

52.4 

1.28 

1.57850 

2.799 

0.564 

54.6 

1.37 

2.12400 

3.389 

0.627 

56.7 

1.46 

2.81050 

4.167 

0.674 

57.3 

1.55 

2.96500 

4.672 

0.635 

50.8 

1.65 

2.39200 

4.964 

0.482 

36.2 

1.75 

1.71825 

5.035 

0.341 

24.2 

1.85 

0.99125 

4.860 

0.204 

13.7 

2.03 

0.05025 

4.977 

0.001 

0.1 

TABLE  E-16 

Photocharacteristics  of  G-8-1  it 23  at  OV  for  H**  Calibration  Data 


X 

I 

H 

R 

V 

(pm) 

(pA) 

(pW> 

(pA/pU) 

(%) 

1.04 

0.00165 

0.1426 

0.012 

1.4 

1.11 

0.01261 

0.9206 

0.014 

1.5 

1.19 

0.03205 

1.9838 

0.016 

1.7 

1.28 

0.04995 

2.8331 

0.018 

1.7 

1.37 

0.06603 

3.3452 

0.020 

1.8 

1.46 

0.08396 

4.0389 

0.021 

1.8 

1.55 

0.08838 

4.5057 

0.020 

1.6 

1.65 

0.07459 

4.7326 

0.016 

1.2 

1.75 

0.05544 

4.8104 

0.012 

0.8 

1.85 

0.03171 

4.6224 

0.007 

0.5 

2.03 

0.00156 

4.7196 
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TABLE  E-17 


Photocharacteristics  of  G-8-1  #23  at  5V  for  H**  Calibration  Data 


k 

I 

H 

R 

V 

Oim) 

(ma) 

(M  W) 

(M  A/MW) 

(%) 

1.04 

0.03550 

0.1426 

0.249 

29.7 

1.11 

0.30600 

0.9206 

0.332 

37.1 

1.19 

0.81250 

1.9838 

0.410 

42.7 

1.28 

1.30950 

2.8331 

0.462 

44.8 

1.37 

1.78400 

3.3452 

0.533 

48.3 

1.46 

2.34100 

4.0389 

0.570 

49.2 

1.55 

2.49250 

4.5057 

0.553 

44.3 

1.65 

2.03050 

4.7326 

0.429 

32.2 

1.75 

1.45900 

4.8104 

0.303 

21.5 

1.85 

0.80850 

4.6224 

0.175 

11.7 

2.03 

0.03550 

4.7196 

0.008 

0.5 

TABLE  E-18 

Photocharacteristics  of  G-8-1  #23  at  10V  for  H**  Calibration  Data 


k 

0m) 

I 

(M) 

H 

OiW) 

R 

(ma/mw) 

V 

(%) 

1.04 

0.04000 

0.1426 

0.281 

33.5 

1.11 

0.31900 

0.9206 

0.347 

38.7 

1.19 

0.85750 

1.9838 

0.432 

45.0 

1.28 

1.37400 

2.8331 

0.485 

47.0 

1.37 

1.84550 

3.3452 

0.552 

50.0 

1.46 

2.44650 

4.0389 

0.606 

51.4 

1.55 

2.58550 

4.5057 

0.574 

46.0 

1.65 

2.11700 

4.7326 

0.447 

33.6 

1.75 

1.48700 

4.8104 

0.309 

21.9 

1.85 

0.83550 

4.6224 

0.177 

11.9 

2.03 

0.03750 

4.7196 

0.008 

0.5 
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TABLE  E-19 


Photocharacteristics  of  D-6-1  #6  at  OV  for  H*  Calibration  Data 


X 

(pm) 

I 

(uA) 

H 

(pw) 

R 

(pA/^W) 

r) 

(%) 

1.04 

0.00083 

0.052 

0.0160 

1.9 

1.11 

0.05565 

0.258 

0.2157 

24.1 

1.19 

0.12995 

0.565 

0.2300 

24.0 

1.28 

0.18515 

0.795 

0.2329 

22.6 

1.37 

0.23055 

0.962 

0.2397 

21.7 

1.46 

0.28080 

1.183 

0.2374 

20.2 

1.55 

0.28205 

1.327 

0.2125 

17.0 

1.65 

0.22895 

1.409 

0.1625 

12.2 

1.75 

0.16650 

1.430 

0.1164 

8.3 

1.85 

0.08895 

1.380 

0.0645 

4.3 

2.03 

0.00360 

1.413 

0.0025 

0.2 

TABLE  E-20 

Photocharacteristics  of  D-7-1  // 6  at  5V  for  H*  Calibration  Data 


X 

I 

H 

R 

R' 

n 

*7 7 

(pm) 

(pA) 

(pW) 

(pA/pW) 

CpX/pU) 

m 

(%) 

1.04 

0.03525 

0.052 

0.6779 

0.409 

81 

57 

1.11 

0.26975 

0.258 

1.046 

0.632 

117 

83 

1.19 

0.71025 

0.565 

1.257 

0.759 

131 

93 

1.28 

1.10850 

0.795 

1.394 

0.842 

135 

99 

1.37 

1.48625 

0.962 

1.545 

0.933 

140 

100 

1.46 

1.95950 

1.183 

1.656 

1.000 

141 

86 

1.55 

2.01025 

1.327 

1.515 

0.915 

121 

57 

1.65 

1.52050 

1.409 

1.079 

0.652 

81 

36 

1.75 

1.02950 

1.430 

0.720 

0.435 

51 

17 

1.85 

0.48400 

1.380 

0.357 

0.212 

24 

5 

2.03 

0.01650 

1.413 

0.012 

0.007 

7 
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TABLE  E-21 


Photocharacteristics  of  D-6-1  #6  at  10V  for  H*  Calibration  Data 


X 

I 

H 

R 

r' 

i 

V 

(finO 

(pA) 

(/iW) 

(ptA/^W) 

(fiA/fiW) 

(%) 

(%) 

1.04 

0.01425 

0.052 

0.274 

0.199 

33 

28 

1.11 

0.22375 

0.258 

0.867 

0.631 

97 

83 

1.19 

0.59125 

0.565 

1.046 

0.761 

109 

93 

1.28 

0.91550 

0.795 

1.152 

0.838 

112 

96 

1.37 

1.22100 

0.962 

1.269 

0.924 

115 

98 

1.46 

1.62500 

1.183 

1.374 

1.000 

117 

100 

1.55 

1.67650 

1.327 

1.263 

0.919 

101 

86 

1.65 

1.26975 

1.409 

0.901 

0.656 

68 

58 

1.75 

0.84500 

1.430 

0.591 

0.430 

42 

36 

1.85 

0.41550 

1.380 

0.301 

0.219 

20 

17 

2.03 

0.01725 

1.413 

0.012 

0.009 

1 

1 

TABLE  E-22 

Photocharacterisitcs  of  D-6-1  #6  at  OV  for  H**  Calibration  Data 


X 

Qiin) 

I 

(l*A) 

H 

(MW) 

R 

(MA/MW) 

r 1 

<%) 

1.04 

0.00755 

0.0405 

0.1864 

14.5 

1.11 

0.04695 

0.2613 

0.1798 

20.1 

1.19 

0.10470 

0.5630 

0.1860 

19.4 

1.28 

0.14660 

0.8041 

0.1823 

17.7 

1.37 

0.18080 

0.9494 

0.1904 

17.2 

1.46 

0.21375 

1.1463 

0.1865 

15.8 

1.55 

0.21145 

1.2788 

0.1653 

13.2 

1.65 

0.17010 

1.3432 

0.1266 

9.5 

1.75 

0.12370 

1.3653 

0.0906 

6.4 

1.85 

0.07260 

1.3119 

0.0553 

3.7 

2.03 

0.00345 

1.3395 

0.0026 

0.2 
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TABLE  E-23 


Photocharacteristics  of  D-6-1  #6  at  5V  for  H**  Calibration  Data 


A 

(.lira) 

I 

(pA) 

H 

(pW) 

R 

(pA/pW) 

R' 

(pA/pW) 

V 

(%) 

v' 

(%) 

1.04 

0.05000 

0.0405 

1.2346 

0.6464 

147 

89 

1.11 

0.35250 

0.2613 

1.3490 

0.7062 

151 

91 

1.19 

0.87800 

0.5630 

1.5600 

0.8167 

163 

98 

1.28 

1.31900 

0.8041 

1.6403 

0.8588 

159 

96 

1.37 

1.73650 

0.9494 

1.8290 

0.9575 

166 

100 

1.46 

2.18950 

1.1463 

1.9101 

1 . 0000 

162 

98 

1.55 

2.21450 

1.2788 

1.7317 

0.9066 

139 

84 

1.65 

1.70150 

1.3432 

1.2668 

0.6632 

95 

57 

1.75 

1.16750 

1.3653 

0.8551 

0.4477 

61 

37 

1.85 

0.62950 

1.3119 

0.4798 

0.2512 

32 

19 

2.03 

0.02600 

1.3395 

0.0194 

0.0102 

1 

1 

TABLE  E-24 

Photocharacteristics  of  D-6-1  #6  at  10V  for  H**  Calibration  Data 


A 

(pm) 

I 

(pA) 

H 

(pW) 

R 

(pA/pW) 

R' 

(pA/pW) 

T] 

(%) 

n‘ 

(%) 

1.04 

0.04450 

0.0405 

1.0988 

0.6182 

131 

85 

1.11 

0.32700 

0.2613 

1.2514 

0.7040 

140 

91 

1.19 

0.81450 

0.5630 

1.4467 

0.8139 

151 

98 

1.28 

1.23100 

0.8041 

1.5309 

,  0.8613 

148 

96 

1.37 

1.61500 

0.9464 

1.7011 

0.9570 

154 

100 

1.46 

2.03750 

1.1463 

1.7775 

1.0000 

151 

90 

1.55 

2.06100 

1.2788 

1.6117 

0.9067 

129 

84 

1.65 

1.O0150 

1.3432 

1.1923 

0.6708 

90 

58 

1.75 

1.08200 

1.3653 

0.7925 

0.4459 

56 

36 

1.85 

0.58500 

1.3119 

0.4459 

0.2509 

30 

19 

2.03 

0.02050 

1.3395 

0.0153 

0.0086 

1 

1 
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TABLE  E-25 


Photocharacteristics  of  D-6-3  #6  at  OV  for  H*  Calibration  Data 


k 

I 

H 

R 

r} 

(turn) 

(H  A) 

<HW) 

(ha/mw) 

(%) 

1.04 

0.01090 

0.052 

0.2096 

25.0 

1.11 

0.06580 

0.258 

0.2550 

28.5 

1.19 

0.14407 

0.565 

0.2550 

26.6 

1.28 

0.20235 

0.795 

0.2545 

24.7 

1.37 

0.24790 

0.962 

0.2577 

23.3 

1.46 

0.29217 

1.183 

0.2470 

21.0 

1.55 

0.28597 

1.327 

0.2155 

17.2 

1.65 

0.22277 

1.409 

0.1581 

11.9 

1.75 

0.15997 

1.430 

0.1119 

6.4 

1.85 

0.08992 

1.380 

0.0652 

4.4 

2.03 

0.00400 

1.413 

0.0028 

0.2 

TABLE  E-26 

Photocharacteristics  of  D-6-3  it 6  at  5V  for  H*  Calibration  Data 


A  I  H  R  R '  T}  r}' 

_ (M)  _  (MW)  (jM/^W)  (/iA//iW)  (%)  (%) 


1.04 

1.11 

1.19 

1.28 

1.37 

1.46 

1.55 

1.65 

1.75 

1.85 

2.03 


0.10000 

0.052 

1.923 

0.64500 

0.258 

2.500 

1.52200 

0.565 

2.694 

2.27300 

0.795 

2.859 

2.94800 

0.962 

3.064 

3.70450 

1.183 

3.131 

3.65950 

1.327 

2.758 

2.79000 

1.409 

1.980 

1.92800 

1.403 

1.374 

0.99800 

1.380 

0.723 

0.04200 

1.413 

0.030 

0.614 

229 

81 

0.798 

279 

99 

0.860 

281 

100 

0.913 

277 

98 

0.979 

277 

98 

1.000 

266 

95 

0.881 

221 

79 

0.632 

149 

53 

0.439 

97 

35 

0.231 

48 

17 

0.010 

2 

1 
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TABLE  E-27 


Photocharacteristics  of  D-6-3  //6  at  10V  for  H*  Calibration  Data 


X 

I 

H 

R 

R’ 

n 

n* 

Qxm) 

(HA) 

(nw) 

(^A/mW) 

QiA/pW) 

(%) 

(%) 

1.04 

0.10350 

0.052 

1.990 

0.629 

237 

83 

1.11 

0.64650 

0.258 

2.506 

0.793 

280 

98 

1.19 

1.54400 

0.565 

2.733 

0.864 

285 

100 

1.28 

2.29950 

0.795 

2.892 

0.915 

280 

98 

1.37 

2.97500 

0.962 

3.093 

0.978 

280 

98 

1.46 

3.74050 

1.183 

3.162 

1.000 

269 

94 

1.55 

3.68800 

1.327 

2.779 

0.879 

222 

78 

1.65 

2.80750 

1.409 

1.993 

0.630 

150 

53 

1.75 

1.94850 

1.430 

1.363 

0.431 

97 

34 

1.85 

1.02050 

1.380 

0.739 

0.234 

50 

18 

2.03 

0.04250 

1.413 

0.030 

0.009 

2 

1 

TABLE  E-28 

Photocharacteristics  of  D-6-3  #6  at  OV  for  H**  Calibration  Data 


X 

(/im) 

I 

O' A) 

H 

(MW) 

R 

QiA/MW) 

V 

(%) 

1.04 

0.00470 

0.0405 

0.1160 

13.8 

1.11 

0.03203 

0.2613 

0.1226 

13.7 

1.19 

0.07328 

0.5630 

0.1302 

13.6 

1.28 

0.10673 

0.8041 

0.1327 

12.9 

1.37 

0.13340 

0.9494 

0.1405 

12.7 

1.46 

0.15908 

1.1463 

0.1388 

11.8 

1.55 

0.15238 

1.2788 

0.1192 

9.5 

1.65 

0.11125 

1.3432 

0.0828 

6.2 

1.75 

0.07573 

1.3653 

0.0555 

3.9 

1.85 

0.04243 

1.3119 

0.0323 

2.2 

2.03 

0.00158 

1.3395 

0.0012 

TABLE  E-29 


Photocharacteristics  of  D-6-3  // 6  at  5V  for  H**  Calibration  Data 


X 

I 

H 

R 

R’ 

>7 

(jim) 

(ma) 

(liW) 

0iA/pW) 

QiA/^iW) 

(%) 

(%) 

1.04 

0.10650 

0.0405. 

2.6296 

0.7106 

313 

96 

1.11 

0.73350 

0.2613 

2.8071 

0.7586 

314 

96 

1.19 

1.76000 

0.5630 

3.1261 

0.8448 

326 

99 

1.28 

2.63250 

0.8041 

3.2738 

0.8847 

317 

97 

1.37 

3.43400 

0.9494 

3.6170 

0.9747 

327 

100 

1.46 

4.24200 

1.1463 

3.7006 

1.0000  ' 

314 

96 

1.55 

4.14500 

1.2788 

3.2413 

0.8759 

259 

79 

1.65 

3.10900 

1.3432 

2.3146 

0.6255 

174 

53 

1.75 

2.14100 

1.3653 

1.5682 

0.4238 

111 

34 

1.85 

1.16550 

1.3119 

0.8884 

0.2401 

60 

18 

2.03 

0.04450 

1.3395 

0.3322 

0.0898 

2 

1 

TABLE  E-30 

Photocharacteristics  of 

D-6-3  #6  at 

10V  for  H** 

Calibration 

Data 

X 

I 

H 

R 

R' 

V 

V' 

(/im) 

(M) 

<^W) 

(ptA/fiV!) 

(/iA/MW) 

(%) 

(%) 

1.04 

0.11050 

0.0405 

2.728 

0.7548 

325 

100 

1.11 

0.72700 

0.2613 

2.782 

0.7698 

311 

96 

1.19 

1.75400 

0.5630 

3.115 

0.8619 

325 

100 

1.28 

2.61400 

0.8041 

3.251 

0.8996 

315 

97 

1.37 

3.40400 

0.9494 

3.585 

0.9920 

325 

100 

1.46 

4.14300 

1.1463 

3.614 

1.0000 

30" 

95 

1.55 

4.22850 

1.2788 

3.307 

0.9157 

265 

82 

1.65 

3.09200 

1.3432 

2.302 

0.6370 

173 

53 

1.75 

2.14350 

1.3653 

1.570 

0.4344 

111 

34 

1.85 

1.15750 

1.3119 

0.882 

0.2441 

59 

18 

2.03 

0.04500 

1.3395 

0.336 

0.0930 

2 

1 

( 
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